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ABSTRACT 



The radioactive decay energy (RDE) deposition in supernovae from the decay chain "'''Ni ^'^^Co ^•''^Fe 
usually directly powers the ultraviolet/optical/infrared (UVOIR) bolometric luminosity of supernovae in 
their quasi-steady state phase until very late times. The result for this phase is exponential/quasi- 
exponential UVOIR bolometric light curves and often exponential/quasi-exponential broad band light curves. 
A presentation is given of a simple, approximate, analytic treatment of RDE deposition that provides 
a straightforward understanding of the exponential/quasi-exponential behavior of the UVOIR bolometric 
luminosity and a partial understanding of the exponential/quasi-exponential behavior of the broad band 
light curves. The treatment reduces to using a normalized deposition function N-^^(t) as an analysis tool. 
(The absolute deposition is determined by specifying the initial ^^Ni mass or fitting absolute supernova 
UVOIR bolometric luminosity.) The time evolution of N^^{t) is determined by three time scales: the half- 
lives of ^^Ni and ^^Co, and a fiducial time parameter to that governs the time-varying 7-ray optical depth 
behavior of a supernova. The to parameter can be extracted from a structural supernova model, or by fitting 
either the RDE deposition curve from a more detailed treatment of deposition or the observed UVOIR 
bolometric light curve from the quasi-steady state phase of a supernova. A to parameter obtained from 
observations can provide a constraint on the important physical parameters of the supernova. The effective 
use of this constraint requires having an adequate parameterized structural supernova model. 

The N^i{t) function is used to analyze the preliminary UVOIR bolometric light curve of SN Ic 1998bw 
(the possible cause of 7-ray burst GRB980425). The SN 1998bw fiducial time to is found to be 134.42 days 
and a prediction is made for the evolution of the SN 1998bw RDE deposition curve and quasi-steady state 
UVOIR bolometric light curve out to day 1000 after the explosion. A crude estimate (perhaps a factor of a few 
too small) of the SN 1998bw mass obtained from a parameterized core-collapse model and to = 134.42 days is 
4.26 Mq. As further examples of the simple analytic treatment, the RDE deposition and luminosity evolution 
of SN la 1992A and SN II 1987A have also been examined. 

The simple analytic treatment of RDE deposition has actually existed for 20 years at least without, 
apparently, being discussed at length. The main value of this paper is the explicit, detailed, general 
presentation of this analytic treatment. 



Subject Headings: supernovae: general — supernovae: individual (SN 1992A, SN 1987A, SN 1998bw) — 
radiative transfer — 7-rays: theory 
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1. INTRODUCTION 



Radioactive beta decay energy (RDE) from radioactive "''^Ni synthesized in the explosion and its daughter 
^^Co is a major source of supernova luminosity. In Type la supernovae (SNe la), RDE powers the whole 
observable luminosity from a few days after explosion (e.g., Harkness 1991) until as far as observations extend 
(i.e., about 944 days after explosion for the case of SN 1992A) or so it seems (Cappellaro ct al. 1997). In 
core-collapse supernovae (i.e., SNe II, lb, Ic), RDE is the main source of luminosity starting from a few days 
to some tens of days after explosion (e.g., Woosley 1988; Young, Baron, & Branch 1995) and perhaps lasting 
until of order 1000 days after explosion (e.g., Simtzcff et al. 1992). Before the RDE becomes dominant, 
shock heat left over from the explosion is the main source. The transition time between shock heat and 
RDE powering depends on the parameters of the supernova, most importantly mass, kinetic energy, density 
profile, and composition. Circumstellar interaction and perhaps pulsar action can also supply energy to 
core-collapse supernovae. 

In this paper we provide a simple, analytic, normalized deposition function to predict RDE deposition 
in supernovae as a function of time. The treatment provides a straightforward understanding of the 
ultraviolet/optical/infrared (UVOIR) exponential/quasi-exponential light curves and a partial understanding 
of the exponential/quasi-exponential broad band light curves of supernovae that occur tens to hundreds 
of days after explosion. The treatment is complementary to the detailed modeling of RDE deposition 
usually performed by Monte Carlo calculations (e.g., Hoflich, Khokhlov, & Miiller 1992; Swartz, Sutherland, 
& Harkness 1995) or, at a less detailed level, with grey (frequency-integrated) radiative transfer (by Monte 
Carlos or by the radiative transfer equation) (e.g., Colgate, Petschek, & Kriese 1980a, b; Swartz et al. 1995; 
Cappellaro et al. 1997; Jeffery 1998a, b). 

The analytic treatment of RDE deposition has previously been presented by Colgate et al. (1980a, b), but 
only for a constant density supernova model and only briefly. We are not aware of any longer presentations. 
Here we present the treatment in detail and in general: i.e., we allow for any sort of homologously expanding 
supernova model. Special cases of supernova models are presented too. 

In § 2 of this paper we briefly discuss the physics of RDE deposition and exponential/quasi-exponential 
light curves. Section 3 develops our RDE deposition function and parameterized structural models for SNe la 
and core-collapse supernovae. As examples of the simple analytic treatment of RDE deposition, we examine 
the RDE deposition and luminosity evolution of SN la 1992A, SN II 1987A, and SN Ic 1998bw in §§ 4, 
5, and 6, respectively. Conclusions appear in § 7. Appendix A presents some useful analytic results for 
homologously expanding supernova models with exponential density proflles. 

2. THE RDE DEPOSITION 
AND EXPONENTIAL/QUASI-EXPONENTIAL LIGHT CURVES 

The RDE comes principally in the form of 7-rays and positrons not counting neutrinos. (Neutrinos 
almost entirely escape the supernova ejecta and do not contribute significantly to the RDE deposition.) The 
7-ray energy that does not escape the siipernova is converted to fast electron kinetic energy almost entirely 
by Compton scattering. The positrons lose nearly all their kinetic energy slowing down by collisions and 
then annihilate to form 7-rays. The lost kinetic energy and the annihilation 7-ray energy (if it does not 
escape) become fast electron kinetic energy also. There is also a small amoimt of RDE in X-rays from inner 
atomic orbital transitions and in the kinetic energy of ejected fast atomic electrons: both forms being a 
direct consequence of the nuclear decay. The X-ray energy (if it does not escape) is also converted to fast 
electron kinetic energy by ionizations. Although the fraction of decay energy in X-rays is very small, it is 
likely to become significant at very late times when Compton optical depth has become very small due to 
decreasing density, but X-ray optical depth is still large because of the large X-ray opacity. The ratio of the 
effective absorption opacities of the decay X-rays to the decay 7-rays (using 0.01 McV as the dividing line 
between X-rays and 7-rays) is of order lO'^-lO"'. The decay-ejected fast atomic electrons lose their kinetic 
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energy to other electrons, of course. In this paper we will lump the positrons and ejected atomic electrons 
together as positron-electron (PE) particles. 

The overwhelmingly dominant decay chain for the observable epoch of most supernovae is 
56Nj ^56cq _^56pg ^^^^ half-lives of 6.077 ± 0.012 days and 77.27 ± 0.03 days for the first and second 
decays, respectively. (Half-life data in this paper is from Lawrence Berkeley National Laboratory Isotopes 
Project Web Data Base 1999, hereafter LBL.) The first decay releases almost all its non-neutrino energy in 
the form of 7-rays and the second in 7-rays and, in 19 % of the decays, in positrons (e.g., Browne & Firestone 
1986; Huo 1992). The ^^Ni and ^^Co 7-ray and mean positron kinetic energy are in the range ~ 0.16-3.6 MeV 
(e.g., Browne & Firestone 1986; Huo 1992). The long-lived radioactive species ^'^Co (half-life 271.79 ± 0.09 
days), ^^Fe (half-life 2.73 ± 0.03 years), and '*^Ti (half-life 63 ± 3 years) are likely to become important at 
later times: hundreds to thousands of days after the explosion. These species (or their short-lived parents) 
are synthesized in comparatively small abundance in the explosion, but their long half-lives gives them 
importance at very late times. 

The RDE converted into fast electron kinetic energy is said to be deposited. The fast electrons caused 
by the deposition ionize, excite, and heat the low-ionization-state supernova plasma. The deposited RDE 
is ultimately converted mainly into UVOIR radiation that escapes the supernova and forms the observed 
supernova light curves. (The complicated cascade process by which the fast electron kinetic energy gets 
converted into other forms of energy is discussed by, e.g., Fransson [1994, p. 688ff] and Liu & Victor [1994].) 
Because radioactive decay is exponential, there is a prima facie reason to expect supernova light curves 
to decline with time at least partially exponentially. Historically, it was the exponential (really quasi- 
exponential) late-time light curves of SNe la that first suggested radioactive sources for these supernovae 
(Baade et al. 1956). 

Supernova light curves cannot be nearly exactly exponential in most cases for a number of reasons. For 
core-collapse supernovae there is the aforesaid early shock heat source, and circumstellar interaction and 
pulsar sources. For all supernovae, the UVOIR diffusion time scale starts out much longer than the dynamic 
time scale and the radioactive decay time scale. Thus, deposited RDE and shock heat energy is initially 
largely trapped and is released increasingly rapidly as the supernova expands and decreases in density. This 
early trapping rules out seeing the signature in the light curves of the ^^Ni decay entirely and the early 
signature of the ^^Co decay. In the case of SNe la, the transition to quasi-exponential light curve decline at 
about 60 days after explosion (or about 40 days after maximum light [i.e., B maximum]) (e.g., Leibundgut 
1988; Leibundgut et al. 1991b) shows that the UVOIR trapping and the time period for processing RDE into 
UVOIR emission have become small. Thus a quasi-steady state has been established (or perhaps slightly 
later at 70 days after explosion [Pinto & Eastman 1996]) in which at any instant the entire radiative transfer 
of the supernova can be treated by a time-independent calculation to good accuracy. In the case of core- 
collapse supernovae, the transition to the quasi-steady state (exponential/quasi-exponential phase of the 
light curves) can occur at various times depending on the nature of progenitor, but for a massive progenitor 
typically at of order a hundred days after explosion as evidenced by SN 1987A (e.g., Bouchet et al. 1991; 
Suntzeff et al. 1991). 

The late-time light curves (i.e., those from after the establishment of the quasi-steady state) are not 
usually truly exponential despite the effectively instantaneous conversion of RDE to UVOIR luminosity for 
two reasons. First, the 7-ray optical depth of the ejecta can cease to be completely trapping before or at 
about the same time the UVOIR optical depth becomes small. The escape of 7-rays from the ejecta adds 
a non-exponential factor to the RDE deposition (as we will show in §§ 3.1 and 3.2). Second, one observes 
light curves in particular frequency bands. Only the late-time UVOIR bolometric light curve in the case 
of complete 7-ray and PE particle trapping (or complete 7-ray escape and complete PE particle trapping) 
is guaranteed to be exponential. The individual bands receive the RDE after a more or less complicated 
time-dependent distribution by atomic processes. 

Nevertheless, the late-time light curves of supernovae, particularly SNe la in the B and V bands, often 
appear very exponential at least in some phases when they are not in fact exponential to 1st order in 
those phases. We define a 1st order exponential phase to one where the logarithmic slope of a function or 
equivalently its instantaneous half-life (or instantaneous e-folding parameter) is constant to first order at 
some point in the phase. We use the expression quasi-exponential to describe phases of fmictions which 
appear very exponential without being 1st order exponential. We will adopt the convention that only light 
curve phases which are 1st order exponential can be considered truly exponential phases. 
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With older light curve data, it was often possible to fit straight lines (on semi-logarithmic plots) within 
the uncertainty to quasi-exponential light curves. For example, Doggett & Branch (1985) fit a straight line 
with a half-life of 44 days to a compilation of SN la blue magnitude data for days ~ 120-320 after explosion: 
the data had a dispersion of order 1 magnitude (0.4 dex) about the line. (Here "blue" refers to a mix of 
bands from various older and newer magnitude systems that sample the blue side of the optical.) Barbon, 
Cappellaro, & Turatto (1984) found a half-life of 50 days for a line fit to another compilation of SN la blue 
magnitude data for days ^ 220-420 after explosion: the magnitudes had a dispersion of order 0.5 magnitude 
(0.2 dex) about the line. Kirshner & Oke (1975) fit a straight line to an AB (roughly B — 0.2) light curve 
for SN la 1972E with a half-life of 58 days for ~ 98-732 days after explosion. (We assume a rise time to 
maximum light of 18 days for SNc la: sec § 4.) Older late-time blue light curve data for SNc la generally 
showed a half-life of ~ 56 days (e.g., Rust, Leventhal, & McCall 1976). This order of half-life originally 
suggested the spontaneous fission of ^^^Cf (half-life 60.5 ± 0.2) as a dominant energy source (Baade et al. 
1956). The ^^Ni ^^^Co ^^^^Fc beta decay chain for SNc la was put forward by Pankey (1962) and Colgate 
& McKee (1969). This RDE source was later accepted for core-collapse supernovae too. 

Very accurate modern data shows that SNe la late-time light curves are not truly exponential, but have 
slowly decreasing logarithmic decline rates (i.e., increasing instantaneous half-lives) (e.g.. Cappellaro et al. 
1997). (That the older SN la decline half- lives cited above increase with later or longer time coverage 
is probably not a coincidence.) The increasing (instantaneous) half-life is easily understood qualitatively. 
The half-life falls below the "''^Co half-life because of 7-ray escape. But as the 7-ray escape increases with 
decreasing optical depth, the fraction of RDE deposition from the more strongly trapped PE particles 
increases and this eventually drives the half-life back toward the ^^Co half-life asymptotically. Colgate et al. 
(1980a, b) quantitatively demonstrated the increasing RDE deposition half-life in numerical calculations. 
More recently Cappellaro et al. (1997) have done the same demonstration with light curve fits to modern 
SN la data. In § 3.2, we show how both decreasing and increasing RDE deposition half-lives can arise 
analytically. In § 3.3, we show why the signature of increasing half-life is what one sees in the late-time 
light curves of SNe la. Both decreasing and increasing half-life phases can (but not necessarily will) have a 
signature in core-collapse supernova light curves (see § 3.4). 

Positrons and very fast electrons can also move in supernova ejecta and can escape. They are in fact 
thought to be kept from traveling far by magnetic fields: possibly tangled magnetic fields. Positron transport 
may be an important process in supernovae and it is being actively investigated (e.g., Colgate et al. 1980a, 
b; Chan & Lingenfelter 1993; Cappellaro et al. 1997; Ruiz-Lapuente 1997; Milne, The, & Leising 1997: Ruiz- 
Lapuente & Spruit 1998). However, it is certain that positrons and fast electrons are much more strongly 
trapped than 7-rays. Since our interest in this paper is in studying the most basic part of RDE deposition 
analytically, we make the simplifying assumption that PE particle kinetic energy is deposited local to the 
creation of the PE particles. 

At very late times, hundreds of days after explosion, the quasi-steady state is predicted to breakdown 
(e.g., Axelrod 1980, p. 48; Fransson & Kozma 1993; Pransson, Houck, & Kozma 1996). This is because the 
time scale for the fast electrons to lose energy and time scale for recombination cease to be short compared 
to the decay and dynamic time scales. Thus deviation from true 1st order exponential decay or quasi-steady 
state quasi-exponential decay in supernova light curves is to be expected for this reason if no other. There 
is also an observational problem at very late times. At some point an "infrared catastrophe" should occur 
(e.g., Axelrod 1980, p. 70; Pransson et al. 1996) in which the bulk of the emission shifts from the optical 
to the infrared where it is usually much less observable. Thus even if the UVOIR bolometric luminosity 
remains exponential/quasi-exponential, the observable light curves might show sudden, sharp declines. 

3. THE RDE DEPOSITION FUNCTION 

In the present case we arc interested only in the global RDE deposition per unit time (i.e., the deposition 
function) for a supernova principally as a means for understanding supernova light curve behavior. First, 
we will present a fairly exact, general expression Dg{t) for deposition and show how it depends on time 
(§ 3.1). Then in § 3.2 we will derive N^^{t): an approximate, analytic, normalized deposition function for 
^^Ni and its daughter ^^Co. It is N-^^{t) and its subcomponent the absorption function f{x) that are useful 
for analytic insight into the RDE deposition and light curve behavior. In § 3.3, we posit a parameterized 
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SN la model and in § 3.4 a parameterized core-collapse supernova model to which N^i{t) is applied. 

3.1. The Function Dg{t) 

A great thing about 7-ray transfer in supernovae is that the non-scattering component of the source 
function (i.e., the radioactive decay source function) is independent of 7-ray transfer. This means that 
detailed 7-ray transfer with Compton scattering can be treated by a non-iterative Monte Carlo calculation. 
But a second great thing is that the actual complicated process of Compton scattering and photon 
degradation in energy can be approximated by a grey, pure absorption opacity treatment with an uncertainty 
of only a few percent at most in local and global energy absorption (Colgate et al. 1980a, b; Sutherland 
& Wheeler 1984; Ambwani & Sutherland 1988; Swartz et al. 1995; Jeffery 1998a, b). This means that 7-ray 
transfer and absorption can be reduced to simple integration and that is what we will do. We also will assume 
that the 7-ray transfer time can be treated as time-independent: i.e., the time scale for 7-ray transfer is 
much shorter than the supernova dynamical time scale and decay time scale. This assumption can begin to 
fail hundreds of days after explosion (e.g., Jeffery 1998b). 

The radiative transfer of X-rays must be treated separately from 7-ray transfer, since X-rays face a much 
higher absorption opacity than 7-rays. But since the small amount of energy in X-rays is only important at 
very late times (hundreds of days after explosion) we will not consider them further in this section: i.e., we 
will not include them in the Dg{t) or N^i{t) functions. 

With aforesaid simplifications an expression for Dg{t) suitable for our purposes can just be written down 
from simple radiative transfer: 

^gW=E / d^rp{f,t)eii^hj, + U^<f^ [l_exp(-T)]| , (1) 

where the sum is over radioactive species i, the first integral is over all volume, p is mass density, a is the 
radioactive species RDE production (not counting neutrinos which are all assumed to escape) per unit time 
per unit mass, pei is the same except per unit volume, /pe is the fraction of RDE (not counting neutrino 
energy) that goes into fast PE kinetic energy (assumed to be all locally deposited), /ph is the fraction of RDE 
(counting positron annihilation energy, but not counting neutrino energy) that goes into photons, the second 
integral averages over all solid angle for each point r, t is the effective 7-ray absorption opacity optical depth 
in a given direction from r to the surface of the supernova (i.e., the optical depth of the beam path), and 
(1 — e~^) is the absorption probability for photons that travel from r to the surface. 

Supernovae after very early times are in homologous expansion where the velocities of all mass elements 
are constants, and the internal gas and bulk kinetic energy lost to P dV work is negligible. (Exceptions to 
homologous expansion probably exist and have important consequences [e.g., Woosley 1988], but we will not 
consider these in this paper.) The radial position r of a mass element in homologous expansion is given by 

f=vt, (2) 

where v is the mass element velocity and t is the time since explosion: the initial radii of the mass elements 
are negligible after very early times. Thus velocity can be used as a comoving coordinate. The density at 
any velocity decreases as t~^. Therefore we can write 

p{f,t) = p{v,t) = po{v) (^j-^ , (3) 

where po{v) is the density at v at a fiducial time to. Note we are not assuming spherical symmetry for the 
supernova. 

Using the homologous condition equations (2) and (3), equation (1) can be rewritten as 

Dsit) = E / '^'^ 4po{v)ei |/PE + /ph / ^ [1 - exp (-r)] | , (4) 
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where the first integral is now over all velocity space volume. The deposition function now can be seen to 
depend on time in only two ways. First, through Ci which in general is a sum of the exponential terms for 
radioactive species i. For beta decay from a species synthesized in the explosion 

e, =Xi'"Ciexp(-Vie,i) , (5) 

where X"" is initial mass fraction of species i, Ci is an energy generation rate coefficient, and te.i is e-folding 
time (i.e., half-life divided by ln[2]). For beta decay from a species whose parent was synthesized in the 
explosion 

= [exp {-t/te,i) - exp (-Vie,pa(i))] , (6) 

where pa(i) identifies the radioactive parent of radioactive species i, X^.^ is initial mass fraction of the 

parent, Bi is an energy generation rate coefficient, te,i is e-folding time for species i, and te,pa(j) is e-folding 
time for the parent. The C and B coefficients are given by 

Q= and B,= ^p^+J^'' ^ (7) 

Tnaniu^ite,i 'Tlamu^pa(i) \te,i ~ te,pa.{i) } 

where Qph+PE.i is the mean photon plus fast PE kinetic energy per decay, Ai is atomic mass, ^pa(i) is the 
parent's atomic mass, and rriamu is the atomic mass unit (amu). For our simplified deposition treatment 
introduced in § 3.2, we consider only ^^Ni and ^^Co decays since these are overwhelming the most important 
until very late times, hundreds of days after explosion. The parameters for ^^Ni and ^^Co decays are given 
in Table 1. 

The second time dependence of Dg{t) comes through the optical depths of the beam paths: 

T= dsKp[v{s),t] = (-^] / dVstoKpo[viVs)] , (8) 

J cm \ / J cm 

where we have used the homologous expansion condition equations (2) and (3) to get the second expression, 
"em" and "sur" specify the point of emission and surface, respectively, in either space or velocity, s and Vg 
are beam path length in physical space and velocity space, respectively, and k is the effective absorption 
opacity. Note the beam path is general: i.e., not usually radial. As can be seen from equation (8), the optical 
depth between any two points in velocity space decreases as t~'^. 

Swartz et al. (1995) showed for the 7-rays important in supernovae that 

K W Kcll~^ . (9) 



The overwhelmingly dominant 7-ray opacity in supernovae is Compton scattering (e.g., Swartz et al. 1995; 
Jeffery 1998b). Thus the opacity is largely determined by the number density of electrons per unit mass and 
this accounts for the factor of the inverse of the mean atomic mass per electron 

Me-^=E^' (10) 



where the sum is over all elements i, Xi is mass fraction of element i, Zi is nuclear change (since all electrons 
are approximately free electrons to the 7-rays important in supernovae), and Ai is again atomic mass. For 
^^Co, the Kc parameter ranges between ^ 0.05 to 0.065 cm^ g^^ and for '"'^Ni, between ^ 0.06 to 0.1 cm^ g^^ 
(Swartz et al. 1995; Jeffery 1998a, b). The variation of the Kc parameter depends on the global optical depth 
structure of the supernova: the low values are for the optically thin limit; the large values for the optically 
thick limit. (Note Kc actually has a very weak composition dependence too that we neglect here.) We will 
not try to incorporate the weak, but complex, global optical depth structure dependence (and hence time 
dependence) of Kc or variations in pe in our formalism directly. We will simply try to choose appropriate 
values of k for the particular cases we examine and treat those values as constants for all locations and 
epochs. 
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Equation (4) woiilci be vastly simplified if we could replace the location-, direction- and time-dependent 
optical depth by a mean optical depth that was only time-dependent. Given all our assumptions about 7-ray 
transfer, an exact mean f would be obtained from 

j d?vpo{v)eij ^ [l-exp(-T)]. 
l-exp(-f) = ^ . (11) 



Evaluating this exact mean f is, of course, tantamount to solving the exact deposition problem. We note 

that f has a complex time dependence through the e^'s and the exponentials of the individual optical depths, 
not the simple t~'^ time dependence of an individual optical depth given by equation (8). In the optically 
thin limit equation (11) reduces to 




(12) 



In this limiting case when there is a only a single radioactive species, f does have the simple t~'^ time 
dependence since the individual Tj's occur linearly in equation (12) and the radioactive decay time dependence 
cancels. 



3.2. The Function N*i{t) 

So far we have developed, given the assumptions we have made, an accurate general expression for RDE 
deposition, Dg{t) (see § 3.1). The normahzed version of Dg{t) restricted to having only ^^Ni at time zero is 
denoted N^i{t). We will now develop an approximate version of AfNi(i), denoted iVf|J;(i), by approximating 
some of the components of equation (4) for Dg{t) (restricted to "''''Ni at time zero) and normalizing. Since 
only a trace of ^^Co is produced in a supernova explosion and other radioactive species are unimportant 
until hundreds of days after the explosion, the assumption that there is only ^^Ni at time zero is excellent 
until very late times. The absolute deposition in a supernova model is, of course, predicted from N^^{t) or 
Nm{t) with the scale set by specifying the initial ^^Ni mass or by fitting the absolute UVOIR luminosity of 
an observed supernova. 

We assume that ''^Ni decays to negligible abundance while the optical depth is still very large: viz., the 
^^Ni 7-rays are completely trapped. This is a fairly good assumption for all supernovae. It begins to fail 
for SNe la which become optically thin rather quickly. But by the time the SN la quasi-steady state phase 
(which is our main interest) has begun, the ^^Ni is negligible. Given the complete trapping assumption for 
^^Ni 7-rays, we replace 1 — exp(— r) for ^^Ni by 1. 

Note that the effective absorption opacity and thus optical depth scale for ^^Ni 7-rays are larger by a 
factor of about 1.5 in the optically thick limit than those for ^^Co (see § 3.1). But in qualitative discussion 
below, we will treat the two optical depth scales as being same. 

For ^^Co, we approximate the time-, location-, and direction-dependent optical depth r given by 
equation (8) by a location- and direction-independent characteristic optical depth 

2 

= x-^ , (13) 

where Tch,o is the characteristic optical depth at the fiducial time to (i.e., is the fiducial [characteristic] optical 
depth) and 



is a reduced time. 



Tch = Tch,0 
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Since the optically thin f has the same time dependence as Tch (assuming only ''^Ni at time z(;ro). we 
make the optically thin f our choice for Tch for SNe la and other low- mass/rapidly-expanding supernovae. 
These supernovae can be well observed in the 7-ray optically thin phase. If we can establish the optically 
thin f exactly, then N^^(t) will agree exactly with N-f^i(t) in the optically thin limit. For massive and/or 
slowly expanding core-collapse supernovae, the fully optically thin phase may be later than observations or 
later than the period in which the simple N^i{t) behavior is unperturbed by RDE deposition from sources 
other than ^^Ni and ^^Co. For such supernovae, Tch should probably be chosen to be most characteristic 
of the supernova epoch one is trying to model. In order to determine the Tch one needs to posit supernova 
models. We do this for SNe la and core-collapse supernovae in §§ 3.3 and 3.4, respectively. 

With the aforesaid approximations and dividing through by its value ai t = (which is a maximum), 
equation (4) reduces to our approximate normalized deposition function 

7V^i(t) = exp {-t/te,m) + G [exp {-t/te,Co) - exp {-t/te,m)] f[x{t)] , (15) 

where G = Sco/Cni = 0.184641 (to more digits than significant for numerical consistency) and 

fix) = /PE + /ph [1 - exp i-x-^)] (16) 

is what we call the absorption function (for ^^Co). The values of /pe and /ph for ^^Co are given in 
Table 1. Note that N^i{t) has only one free parameter to\/7"ch,o which relates reduced and real time through 
equation (14). This parameter reduces to just to because we will always set Tch.o = 1 since this conveniently 
makes to roughly the time of transition between the optically thick and thin epochs. 

Figure la shows a logarithmic plot with the deposition function N^i{t) for a range of to values (with 
Tch.o = 1 of course). Figure lb shows a logarithmic plot of the absorption function f{x) and three exponential 
fits (which are linear on a logarithmic plot) that we discuss below. 

In Figure la, the complete 7-ray trapping deposition curve is the curve N^i{t) approaches as to — > 00. 
After the ^^Ni contribution becomes negligible the complete trapping curve decays with the ^^Co half-life. 
The complete 7-ray escape deposition curve is the curve that N^i{t) approaches as to 0. The complete 
escape curve is only a complete escape curve for ^^Co 7-rays: it is a complete trapping curve for ^^Ni 7-rays. 
The late-time asymptotic limit of the complete escape curve and all the curves, except for the complete 
trapping curve, is just G/pe exp (— t/te,Co) which also decays with the ^^Co half-life. Note that the smaller 
to, the more rapidly the ejecta becomes optically thin and faster the N^^{t) approaches the complete escape 
curve. We display the complete trapping and escape curves on all subsequent deposition figures for convenient 
reference. 

The N^i{t) curves for to<te,Ni days are not physically realistic. Supernovae with such small to values 
(and there is no evidence that there are such supernovae) would become optically thin before the ^^Ni 
contribution had become negligible. We have assumed that ^^Ni becomes negligible before the complete 
trapping phase is over. We merely display the small to cases to help demonstrate the trend with varying to- 

An exact exponential cannot be a sum of two terms, unless they have the same e-folding parameter. 
Thus iV^j(i) cannot be an exact exponential. However, it very closely approximates an exponential in some 
cases. To dispense with it at once, consider the unphysical case (given our assumptions) where to ;Sie,Ni- In 
this case, the second term of equation (15) decreases quickly (after a very short growth phase from zero) and 
for awhile is small relative to the first term (a ^®Ni exponential). Thus, there will period of nearly exponential 
decay with almost the ^^Ni half- life before the ^^Ni exponential rather abruptly becomes negligible and N^j{t) 
closely approaches its late-time asymptotic limit. 

The ■''^'Ni nearly-exponential phase of N^-^{t) is actually not a 1st order exponential phase (see definition 
in § 2) although it can be very exponential for small enough to as Figure la shows. A sum of exponentials 
with different e-folding parameters can only be considered a 1st order exponential if all the exponential terms, 
except the dominant exponential term, are exponentially small compared to the dominant term. By being 
exponentially small, we mean that an exponential term can be approximated by its value at infinity where 
it is a constant zero to all orders in a series expansion about infinity. The non-^®Ni exponential terms in the 
^^Ni nearly-exponential phase are never really exponentially small compared to the ^^Ni exponential term 
even for to 0. In particular, it is easy to show that the instantaneous e-folding time of the complete escape 
curve (the to ^ case of N^i{t)) is never a first order constant until the ^^Ni exponentials are exponentially 
small. 
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Given that "''''Ni is exponentially small, the sufficient condition for N^^(t) being a 1st order exponential 
is that /(x) be either a 1st order constant or a 1st order exponential itself. The f{x) function has two 
stationary points: a maximum of 1 at a; = and a minimum of /pe at a; = c». These give what we will call 
the early and late 1st order exponential phases, respectively. Both, of course, have the half-life of ^^Co. 

In order to realize the early 1st order exponential phase the complete trapping phase for ^^Co 7-rays 
must continue after the ^^Ni contribution is exponentially small: i.e., f{x) must not differ significantly from 1 
until after the ^^Ni contribution is exponentially small. (Note |log[/(a;)]| < 0.001 for x < 0.4 and f{x) > 0.99 
for X < 0.46.) In Figure la, only the N^i{t) deposition curve for to = 200 days has a noticeable early 1st 
order exponential phase. An early 1st order exponential phase may or may not have a direct signature in 
light curves depending on whether or not it lasts into the quasi-steady state period. 

After the ^^Ni contribution has become exponentially small or after the early 1st order exponential 
phase if it exists, N^^{t) will decline more rapidly than an exponential with ^^Co half- life until very late 
times (on the to time scale) (as Fig. la shows) because of increasing 7-ray escape from the ejecta. But at 
very late times the 7-ray escape probability approaches 1 and 1 — exp (— a:~^) x~'^ 0. Thus f{x) will 
asymptotically approach its minimum value at infinity, /pe = 0.0320. In the asymptotic limit, deposition 
will be at "''^'Co decay rate and entirely due to the decay PE particle kinetic energy which we assume to be 
entirely locally deposited: see § 2 for a caveat about this assumption. (Recall also that we are neglecting 
X-rays.) As the asymptotic limit is approached one has a late 1st order exponential phase as the curves with 
smaller to values in Figure la suggest. The approach to the late 1st order exponential phase is slow on the 
to time scale because the fractional contribution of 7-rays at late times declines only as {to/t)^. On the time 
scale of the ^^Co half-life the approach can be slow or fast, of course, as Figure la shows. 

Although there is no sharp transition to the late 1st order exponential phase, the time when the 7-ray 
and PE deposition are equal can be taken as a conventional transition xtr- When the ratio of PE to 7-ray 
deposition is R, the reduced time (assuming the ^^Ni contribution to deposition is negligible) is 



For i? = 1, we have Xtr « 5.454099. At xn, f{x) is 0.301030 dex (or a factor of 2) above its asymptotic value 
and at 2xtr, 0.098007 dex (or a factor of ~ 1.25). 

The slow approach to the late 1st order exponential phase gives rise to a continuum of quasi-exponential 
phases as Figure la suggests. As an example (chosen for reason: see § 4), we have fitted an exponential 
i^cocf <2xp (— x/xe) to /(x) in the x-range [2,3]. The fit is shown in Figure lb and its parameters are given 
in Table 2. As can be seen in Figure lb the fit closely matches /(x) over the x-range [2,3]: the discrepancy 
between the two is always less than 0.01 dex or about 2.5 % in the x-range [2, 3]. 

Would discrepancies from exponential behavior arising from the curvature of the In [/(x)] function be 
seen in late-time supernova light curves? Late-time light curve data for the various bands from before circa 
1980 often had uncertainties of order 0.5 magnitudes or 0.2 dex. Clearly, discrepancies of ~ 0.01 dex or 
even ^0.1 dex between late-time light curves and exponential fits would be hard to notice or assess for 
significance. As one observed for longer, larger discrepancies could arise. But the later the light curve, the 
fainter the supernova, the more the uncertainty in the observation. Modern, accurate late-time light curves 
for SNe la do show discrepancies from exponential behavior that are a least partly due to the non-exponential 
nature of /(x) in its slow approach to its late asymptotic limit (e.g., Cappellaro et al. 1997; see also § 4). 
The logarithmic decline rates of the light curves become smaller with time. It is now recognized that the 
behavior of late-time SN la light curves from say day 60 after the explosion until very late times is only 
quasi-exponential. 

Besides the early and late 1st order exponential phases, N^^{t) predicts that there can be a third 1st 
order exponential phase. This is when /(x) itself becomes 1st order exponential. In Figure lb, it can be 
seen that ln[/(x)] has an inflection point near x = 1 where the logarithmic slope becomes 1st order constant 
and /(x) becomes a 1st order exponential. Assuming the ^®Ni contribution has become exponentially small, 
then before the inflection point the RDE deposition logarithmic decline rate increases (i.e., the instantaneous 
half-life decreases) and after, the RDE deposition logarithmic decline rate decreases (i.e., the instantaneous 
half-life increases). 




(17) 
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The inflection point XinR can be solved for numerically from an iteration relation obtained from the 
second derivative of ln[/(x)]: 



Xi — 



/ph cxp {-xj^ 



(18) 



where i is iteration number. One obtains Xinfl « 1.040765 and / (xinfl) ~ 0.615465. The parameters for an 
exponential fitted to f{x) at the inflection point are given in Table 2. In the x-range [0.8, 1.35], the maximum 
deviation of f{x) from the inflection point fit is ^ 0.004 dcx (0.01 magnitudes). In the x-range [0.56, 1.82], 
however, the maximum deviation has grown to ~ 0.04 dex (0.1 magnitudes). 

An inflection point 1st order exponential phase will be realized in RDE deposition if ^^Ni abundance 
has become exponentially small by (for the sake of dcfinitcness) x « 1.35 (and, of coursc\ provided that 
-^Ni(i) adequately describes deposition). In Figure la, the N^-^{t) deposition curves for values 200, 100, 
and 60 days all show inflection point 1st order exponential phases. The io = 30 day curve does not obviously 
show one. The ■'^'^Ni fraction of RDE production (not deposition) falls below 50 % by day 18, 20 % by day 30, 
7.6 % by day 40 (x = 4/3 for to = 30 day), and 1 % of the RDE at day 60. From these numbers, it is clear 
we would not expect much of an inflection point 1st order exponential phase for to ^ 30 days. 

A direct signature of a 1st order exponential RDE deposition in the light curves appears only if the 
supernova is in quasi-steady state phase. If there is an inflection point 1st order exponential phase entirely in 
the quasi-steady state phase, observers measuring a UVOIR bolometric light curve or a tracer of the same in 
a largish region about the inflection point would probably infer a slower exponential decline than that given 
by the exponential fit at the inflection point. This is because the finite uncertainty in their data would hide 
some of the logarithmic absorption function curvature. In Table 2 we give narrow (n), middling (m), and 
broad (b) inflection region fits to f{x). The maximum discrepancy of these fits from f{x) are ~ 0.014 dex 
(0.035 magnitudes), ^ 0.028 dex (0.07 magnitudes), and ^ 0.04 dex (0.1 magnitudes). All but the very best 
modern supernova photometry might be unable to show systematic deviations of 0.035 or 0.07 magnitudes 
from a straight fit to light curves. Thus a UVOIR bolometric light curve or tracer would probably appear 
very exponential for a period significantly longer than that of the true 1st order exponential behavior. 

The parameters for the exponential fits to f{x) that we have found of particular interest in this paper are 
collected in Table 2. In addition we have included in Table 2 the half-lives of RDE deposition from strictly 
^^Co decay that result from the product of the ^^Co decay exponential with the f{x) function exponential 
fits for the two fiducial supernova models we present in §§ 3.3 and 3.4. The models are needed to convert 
the half-life of a fitted exponential into a real time half-life tf*2- The half- life of a model's RDE deposition, 
t™i'2-i is then given by 

'■1/2 + '■l/2,Co 

where ti/2.co is the "''^Co half-life. 

Finally, we should take stock of the accuracy of N-^^{t) in comparison to N-^i{t). First, we note that our 
characteristic optical depth with its simple time dependence cannot in general be the real mean optical of 
the ejecta discussed in § 3.1. Thus in general N^^{t) cannot be exactly A^Ni(i) at all times. For SNe la and 
other low-niass/rapidly-expanding supernovae, we have already chosen Tch to be the optically thin r. Insofar 
as the optically thin f can be determined exactly, N^-^ (t) with this choice of Tch will be exact in the optically 
thin limit. On the other hand, in the optically thick limit at very early times, N-^^{t) and A'ni should be 
nearly exactly the same for any Tch of order of the true mean optical depth since virtually all 7-rays are 
locally trapped. 

One can plausibly expect that if iVJ;(t) is exactly correct in two opposite time limits, then it interpolates 
to good accuracy in the intervening transition epoch when Tch is of order 1. Nevertheless, iVJ;(t) can be 
expected to be at its worst in this period. The reason is that the real individual optical depths of a supernova 
will probably range over values larger and smaller than 1 giving rise to complex location and direction varying 
7-ray escape and trapping behavior that a simple optical depth prescription cannot easily mimic exactly. 
Without having the real f for the transition epoch (with its in general complex time dependence), deviations 
of A^Nj(t) from iVNi(i) will occur. Of course, a 1st order exponential phase corresponding to the single 
infiection point one predicted by N-^-^{t) can still exist. If the supernova ejecta is smoothly varying in its 
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properties, a smooth transition in RDE deposition from early behavior to a late 1st order exponential phase 
(with ^^Co half-life of course) can yield a single inflection point 1st order exponential phase at some point 
provided the ^^Ni contribution becomes exponentially small soon enough. But if the ejecta has several 
high density clumps or shells, then the RDE deposition curve might show complex behavior with multiple 
inflection point exponential phases that occur as these clumps or shells make the transition from optically 
thick to thin at different times. 

There is a case in which N^i{t) would always be exact if the effective 7-ray absorption opacity n were a 
true constant. In the limit that the ^^Ni is confined to the exact center of a spherically symmetric supernova, 
f is the exactly the radial optical depth to the center at all epochs and hence f always has the simple t~'^ 
time dependence. The choice of the radial optical depth from the center for Tch would then cause N^^ (t) to 
be exactly Nm{t)- Core-collapse supernovae are expected to have their ^^Ni confined to the central region 
unless there is very extensive mixing. Thus, the exact case of N^i{t) could actually be approached if k were 
constant. Unfortunately, dependence of k on the optical depth structure of supernovae and hence on time, 
albeit weak, is not negligible. 

3.3. A Parameterized SN la Model 

In order to determine the optically thin f (our choice for Tch for SNe la and other low- mass/rapidly- 
expanding supernovae) and relate real time t and reduced time x, we need a structural supernova model. In 
this section we will specify a simple parameterized structural model for SNe la. 

Spherically symmetric hydrodynamic calculations of SN la explosions often (but not always) produce 
models with density profiles that are very exponential (i.e., inverse exponential) with velocity after 
homologous expansion has set in. For example, the well regarded Chandrasekhar mass SN la models 
W7 (Nomoto, Thielemann, & Yokoi 1984; Thielemann, Nomoto, & Yokoi 1986), DD4 (Woosley & Weaver 
1994), and M36 (Hoflich 1995, Fig. 10, but note that the density is mislabeled as energy deposition) are 
quite exponential with equivalent-exponential model e-folding velocities (see Appendix A, eq. [AlO]) of 
about 2700 km s"-'^, 2750 km s"-'^, and 3000 km s~^, respectively. Such nearly-exponential density profile 
models have been quite successful in reproducing SN la spectra (e.g., Jeffery et al. 1992; Kirshner et al. 
1993; Hoflich 1995; Nugent et al. 1995). Therefore we will assume a spherically-symmetric, exponential 
density profile model (i.e.. an exponential model for brevity) for our homologous epoch, parameterized SN la 
model. In Appendix A, we present a number of useful analytic results for exponential models and give a 
prescription for exactly exponential models (equivalent-exponential models) that can approximately replace 
nearly exponential hydrodynamic explosion models. 

The density profile of an exponential model (for the homologous epoch) is given by 

p{v,t) = pce,n(^^ exp (-W/Ue) = Pce,0 ^ CXp (-2:) , (20) 

where Pce,o is the central density at fiducial time to, v is the radial velocity, Vg is the e-folding velocity, and z 
is radial velocity or radial position in velocity space in units of Vg- Substituting for density from equation (20) 
into equation (8) and assuming the opacity k is constant, we find for an exponential model that the 7-ray 
optical depth from an emission point z to the surface (which is at infinity) is 

(^^\ 2 „oo 
-j] J dzs exp {-z') , (21) 

where Tce,o is the radial optical depth to the center at the fiducial time to (see eq. [A16] in Appendix A for 
the expression), Zg is beam path velocity length in units of the e-folding velocity, and 



z' = ^/z^+zf+2zzsil . (22) 

The fx is the cosine of the angle at the emission point between the outward radial direction and the beam 
propagation direction. For a beam in the outward radial direction, 1^ = 1 and the optical depth expression 
reduces to ^ 

Tr = rce,ofyj exp (-2;) . (23) 
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Given that there is only ''^Co (which was initially ^^Ni), equation (12) for the optically thin r changes 
for our exponential SN la model to 

t = ree,o(y) ?, (24) 

where the q parameter is defined by 



^ dzz^exp{-z)X'^i J rf^, exp(-^') 



^ J,°°dzzH.p{-z)X^i ■ ^^^^ 

Note that the initial ^®Ni fraction is location-dependent although time-independent, of course. 

A bit of analysis of equation (25) shows that when all ^^Ni is concentrated in the center, q has its 
maximum value of 1 and the optically thin f, its maximum value of Tcc.o (to/t)'^ ■ The q parameter, in 
fact, acts as measure of ^^Ni concentration. The more the ^^Ni is spread to low optical depth regions, the 
more important the small optical depth contributions are in equation (25), and the smaller q becomes. We 
generalize the q parameter beyond the definition of equation (25) (which is specific to exponential models) 
and consider g as a general initial ^^Ni concentration parameter regardless of the exact supernova model. 
The general q is large for high concentration, small for low concentration. 

To study the behavior of the (exponential model) q parameter consider the special case that the (initial) 
^^Ni exists only between velocity radii a and b (in units of the e-folding velocity of course) and that it is a 
constant there. Then equation (25) reduces to 



/ dzz exp (—2;) / / dzsex.p(—z 

Ja J-i 2 Jo 



dz z'^ exp (— z) 



(26) 



We have been unable to analytically solve the triple integral in the numerator of equation (26) or even find 
an accurate analytic approximation: crude approximate analytic solutions can be found for certain special 
cases. The integrations in the numerator can, however, be easily done numerically. The denominator can be 
solved analytically using equation (Al) in Appendix A. 

Equation (26) is actually heuristically useful for SNe la since hydrodynamic models often show the ^^Ni 
is concentrated in a distinct layer where its abundance is overwhelmingly dominant. For example, model W7 
has an initial "''^Ni mass fraction of over 0.5 in the range ~ 3375-9750 km s~^. In most of this range the mass 
fraction is ^ 0.9. Outside of the range the mass fraction falls fairly quickly to very low values. Using the 
e-folding velocity 2700 km s^^ which is in round numbers the equivalent-exponential model e-folding velocity 
for model W7 (see Appendix A, eq. [AlO]), we obtain from equation (26) by numerical integration q = 0.3591. 
The q value extracted from a full 7-ray transfer calculation for the model W7 (assuming model W7 is an 
exponential model with certain fiducial parameter values: see the description below) yields q = 0.3308. The 
difference is caused by the fact that model W7 is not exactly exponential in density (e.g., it has local density 
maxima at ~ 9200 kms~^, ~ 13300 kms~^, ~ 15000 kms~^, and cuts off sharply at ~ 22000 kms~^) and 
the ^^Ni mass fraction is not exactly constant inside of a sharply defined layer. 

If we set a and b in equation (26) to, respectively, zero and infinity (i.e., an evenly-spread-'''^Ni case), 
we obtain a, q = 0.33333326 at the highest accuracy we have computed. This value is so close to 1/3 that it 
is likely that the actual q value for the evenly-spread-^^Ni case is exactly 1/3. We have not, however, been 
able to find an analytic proof. The fact that the q value for the evenly-spread-^^Ni case is not so different 
from the model W7 and approximated model W7 q values cited just above shows that the interior ^^Ni at 
large optical depth dominates the deposition in this case. Because the 1/3 value is a simple rational number 
that seems to be exactly correct for a particularly well defined system and is not so different from values 
likely to be obtained for more realistic SN la models, we adopt 1/3 below as our fiducial q value for our 
parameterized SN la model. 

From equation (24) evaluated at the fiducial time to with Tcc.o substituted for from equation (AI6) in 
Appendix A, we obtain an expression for a fiducial time in terms of a fiducial characteristic optical depth 
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Tch,o (chosen to be the optically thin r at the fiducial time): 

_ Im nq 1 

Y StT Tch,0 Ve 

We will now rewrite equation (27) in terms of fiducial values. At present Chandrasekhar mass C-0 white 
dwarfs are favored SN la progenitors, and so the fiducial mass is chosen to be 1.38 M©. Since we want our 
expression to be as exactly valid as possible in the optically thin phase, wc set the fiducial k to 0.025 cm^ g~^: 
a general good opaeity value for all-metal ejecta {ne ~ 2) in the optieally thin limit (see § 3.1). We take 
Ve — 2700 km s^^ which in round numbers is the equivalent-exponential model e-folding velocity for model W7 
(see Appendix A, eq. [AlO]): model W7 is representative of currently favored SN la models. The fiducial 
q is set to 1/3 for reasons given above. As explained in § 3.2, the fiducial Tch,o is set to 1 so that to is 
eonvenicnitly roughly the time of transition between the optically thick and thin epochs. With these choices 
the expression for fiducial time becomes 

to = 40.895902 days x J ( ( WT^^,^^^-^) f— ) 

^ Y \1.38MqJ V 0.025 cm2g-i; V^ch.o; \l/^J 

^ ^2700kms-i^ 

For numerical consistency here and throughout this paper we treat fiducial values and the solar mass unit 
Mg = 1.9891 X 10^^ g (Lide & Erederikse 1994, p. 14-2) as nearly exact numbers. Our fiducial SN la model 
is the one that has the fiducial values and thus has fiducial time to = 40.895902 days. 

To test our N^i{t) expression and our parameterized SN la model, we have fitted N^^{t) to a normalized 
W7 RDE deposition curve calculated using a grey radiative transfer procedure and the actual model W7. The 
LS grey 7-ray transfer procedure (Jeffery 1998a, b) was used for the calculation. This is a free-parameter- 
free grey procedure that accounts crudely for the variation of effective absorption opacity with optical depth 
structure. The LS procedure has an uncertainty of order 10% at most for total RDE deposition: the 
uncertainty vanishes in the optically thick and thin limits. For a fair comparison to N^^{t), we included only 
the ^®Ni and ^®Co radioactive sources and turned off X-rays in the LS procedure calculation. Because of 
these limitations on the RDE deposition, we call the curve computed by the LS procedure, the limited W7 
deposition curve. 

The fitting was done by choosing a to value (with Tch,o = 1 as always) that made the N-^-^{t) equal 
the limited W7 deposition curve exactly on day 2000 after explosion. The day 2000 epoch is well into the 
optically thin phase, and so the fitted to reproduces the time-varying optically thin f to within numerical 
accuracy. The fitted fiducial time to = 40.74 days. This model W7 fiducial time is very close to the fiducial 
SN la model fiducial time 40.895902 days. Assuming model W7 is exactly an exponential model with the 
prescribed fiducial parameters, except for the to and q values, we invert equation (28) using to = 40.74 days 
to obtain q = 0.3308. This q cannot have exactly the interpretation dictated by the q equations (25) and (26) 
since model W7 is not exactly an exponential model. However, in the general, if vague, interpretation of 
g as a concentration parameter, the q from the fit has a useful meaning. Model W7 has roughly the same 
order of concentration of ^®Ni as the exact exponential models we considered above. 

Figures 2a and 2b display the limited W7 deposition curve and the fitted N^^{t) curve to day 100 and 
day 500 after the explosion, respectively. The N^i{t) curve is in very close agreement with the limited W7 
deposition eurve and consequently is difficult to see on the Figiire 2b scale. The agreement is closest at early 
and late times as expected: less than ~ 0.001 dex discrepancy before about day 10 and after about day 240. 
The only period when the discrepancy is >, 0.01 dex the day 15 73 period where the N^^{t) curve rises above 
limited W7 deposition curve with the maximum disrepancy being ^ 0.07 dex (about 18 %) at day 33. It was 
expected that the largest discrepancies would occur during the transition from the optically thin to thick 
epochs: i.e., during the phase centered on reduced time x — I (see § 3.2). For model W7 this is the period 
centered on about day 40 after the explosion. (Recall the uncertainty in the limited W7 deposition curve is 
also largest in transition period, but it is no more than of order 10 % percent.) 

We have also plotted on Figures 2a and 2b the complete trapping and escape curves discussed in § 3.2. 
As in Figure la, the complete trapping curve shows the transition between the early phase where declining 



(27) 



(28) 
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■''^Ni and increasing "''^Co provide the RDE and the later phase where dec;reasing "''^Co provides the RDE. 
As we discussed in § 3.2, ^^Ni contribution is truly exponentially small only by of order day 60 which is 
~ 1.5to for both the fiducial SN la model and model W7. Thus, these models cannot really have either an 
early-time nor an inflection point 1st order exponential phase in RDE deposition. For an inflection point 1st 
order exponential phase, the ^^Ni contribution must be exponentially small by of order ~ 1.35fo at about the 
latest. Figure 2a, examined closely, confirms the lack of a first order exponential phase centered on about 
to = 40 days in the model deposition curves. Since the fiducial SN la model and model W7 are reasonable 
models for SNe la, we conclude that SNe la are very unlikely to have either an early time or an inflection 
point 1st order exponential phase. 

Recall that the quasi-steady state of SNe la does not start until about day 60 after the explosion. Thus 

even if the early or inflection point 1st order exponential phases in RDE deposition existed, they would not 
have a direct signature in the light curves. Because SNe la are likely already in the post-inflection point epoch 
when the quasi-steady state phase begins, one expects to see the instantaneous half-lives of their light curves 
only increase with time (i.e., the logarithmic decline rates decrease with time) in the quasi-steady state phase 
at least until very late times. And this is in fact what observations suggest happens (e.g., Cappellaro et al. 
1997; see also § 4). Comparing the deposition curves to the complete escape curve in Figure 2b shows how 
gradual the approach is to late 1st order exponential phase on the time scale of observations which have very 
rarely gone beyond 500 days. 

The real time half-lives of the exponential fits to RDE deposition that result from the product of pure 
^^Co exponential decay and the f{x) function exponential fits discussed in § 3.2 and put on the real time 
scale using the fiducial SN la model to are shown in Table 2. As we have seen above, the early-time and 
inflection point 1st order exponential phases in RDE deposition cannot be realized for the fiducial SN la 
model nor model W7 nor in all likelihood actual SNe la. Thus the exponential flts that correspond to those 
cases and similarly the fit to the x-range [0.23,0.53] cannot be realized. The half-lives corresponding to 
those cases are presented pro forma only. However, the exponential fits to RDE deposition for the x-ranges 
[2,3] and [4,8], and for the late-time limit (i.e., to the time when x oo) arc good fits (after choosing the 
appropriate constant coefficient for the exponential fit) to the N^^it) and limited W7 deposition curves since 
they occur after the ^^Ni is exponentially small. The a;-range [2,3] and [4,8] fits correspond to the periods 
^ 80- 120 and ^ 160 320 days, respectively, for both the fiducial SN la model and model W7. But these are 
only exponential fits to quasi-exponential phases. 

Finally, we note we just needed model W7 itself not the parameterized SN la model in order to 
determine to by fitting. The parameterized model did, however, predict model W7's to value fairly well 
from reasonably chosen fiducial values. Furthermore, the parameterized model gives insight into how the 
deposition is determined. 

3.4. A Parameterized Core-Collapse Supernova Model 

Most supernovae la are fairly homogeneous in their behavior. Core-collapse supernovae are much more 
heterogeneous. Thus it is possible that the parameterized exponential SN la model presented in § 3.3 may be 
fairly adequate for most or all SNe la. But a parameterized core-collapse supernova model is likely to apply 
only sometimes. Here we will develop a parameterized model based on the spherically symmetric 16.2 Mq 
explosion model lOH of Woosley (1988): about 14.5 Mq of the model are ejected and about 1.7 Mq are left 
in a compact remnant, presumably a neutron star. Model lOH is able to account for many of the features of 
SN n 1987A. 

The velocity profile of model lOH in the homologous epoch is very roughly dcscribablc as linear with 
respect to interior mass over most of the ejecta mass: the most linear region extends over the range ~ 4 Mq— 
~ 15 Mq corresponding to ~ 1000-4250 km s~^ (Woosley 1988, Fig. 12a). The rough linear relation between 
velocity and mass implies that density goes roughly as the inverse square of velocity over most of the ejecta 
mass. In the innermost region below ~ 3M0 (where most of the mass was left in the compact remnant), 
the model lOH ejecta density decreases inward roughly speaking (Woosley 1988, Fig. 28): there is no density 
singularity at the center of the ejecta. 

For our parameterized model in the homologous epoch, we assume all the mass is located between 
velocities iia and Vb, and that density goes exactly an inverse-square of velocity in this interval. Thus total 
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mass is given by 

3 



M = 47r ^ " dt; vH^'po (^^^ (^) ' = ^Trpo^o^^o (^b - v,) , (29) 

where po is a fiducial density at a fiducial radial velocity vq and fiducial time to. The radial 7-ray optical 
depth from a velocity v to the surface is given by 



t{v) 



f dvtnpo (f (^)^ = t-^npotlvl (i - 1) , (30) 



where we have assumed k is a constant. 

We now parameterize the characteristic optical depth of the ejecta Teh by the expression 

Tch = T{va/q) , (31) 

where q is again a concentration parameter. The larger q, the smaller v^/q, the larger Tch, and the more 
concentrated the initial ^^Ni must be in the central region of the cjccta. The chosen parameterization of 
Tch is, of course, based on the idea that actual mean optical depth will depend on the overall supernova 
parameters in roughly the same way as a radial optical depth from some velocity to the surface. 

We note that if all the initial ^^Ni is concentrated just at the innermost physical layer (i.e., at Wa), the 
mean optical depth for 7-ray deposition is larger than r (va). A beam of 7-rays from a point at Va, must cross 
a larger optical depth than t (va), except in the outward and inward radial directions. Thus q can be larger 
than 1. This contrasts with our parameterized SN la model where 1 was the maximum value of q (see § 3.3). 

Taking r {va/q) for Tch and using equations (29) and (30), we can now solve for the time when Tch is the 
characteristic optical depth: 

/m k [1 - (V^;b) (1/g)] q ,32^ 

y 47rTch [1 - {Va/V\,)] VaVh ' 

We will rewrite equation (32) in terms of fiducial values guided mainly by model lOH and SN II 1987A. 
For mass we use 15 Mq which is roughly the model lOH ejected mass. The most linear region of model lOH 
cited above contains most, but not almost all, the model lOH mass. Thus it is difficult to select Va and 
Vh based model lOH alone without considerable arbitrariness. But to fit the quasi-steady state UVOIR 
bolometric light curve of SN 1987A (see § 5), we select the values Va = 700kms~^ and fb = 5000 kms~^ 
which give a good fit (along with our other fiducial values) while being consistent with what one might 
estimate from model lOH (Woosley 1988, Fig. 12a). The fit is not uniquely good, of course. In model lOH, 
the ^®Ni is confined to the lowest ejecta velocities (Woosley 1988, Fig. 4). This suggests q to be of order 1. 
It is clear, however, that some ^^Ni was mixed outward perhaps as far as halfway or more in mass fraction 
through the cjccta (e.g., Nomoto et al. 1994a, p. 528ff, and references therein). This suggests a q parameter 
less than 1, maybe much less. But for the sake of definiteness we take the fiducial q = 1. We again choose 
Tchfi = 1 so that our fiducial time to is conveniently roughly the time of transition between the optically 
thick and thin epochs (sec § 3.2). 

Model lOH has of order 10 Mq of metals and helium and of order 5 Mq of hydrogen-rich material which 
is presumably solar or sub-solar in metal (Woosley 1988, Figs. 4 and 17). For a mixed composition of this 
matter using Table 1 of Jeffery (1998a) or Table 2 of Jeffery (1998b), we find k « 0.04cm2g-i for the 
optically thick limit and k « 0.03 cm^g^^ for the optically thin limit. The transition, broadly speaking, 
from optically thick to optically thin phase of SN 1987A spans whole of the SN 1987A quasi-steady state 
UVOIR bolometric light curve (see § 5). Thus we take average value of the limiting opacities 0.035 cm^g~^ 
as our fiducial n. 

With the specified fiducial values, the expression for the fiducial time is 
to = 563.969041 days 



M 



15MqJ V0.035cm2g-iy y^^^^^ 



1 - jva/vb) (1/g) 

1 - (Va/Vh) 



700kms-i\ / 5000 km s 



-1 



if) ■ 
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The parameterized rnodc!l with all the fiducial values chosen is our fiducial core-collapse model and 
563.969041 days is its fiducial time. 

In Figure 3 we show the fiducial model N^^ (t) deposition curve. The to value evaluated for the fiducial 
mo(k;l is so large; that N^^{t) will realize all the 1st order exponential and quasi-exponential RDE deposition 
phases discussed in § 3.2, except, of course, for the ^^Ni nearly-exponential phase which requires small to- 
This is because ^^Ni contribution to the RDE deposition will become exponentially small by of order day 60 
after explosion (as discussed in § 3.2 and as best seen in Figure 2a), and at this time the ejecta is practically 
speaking still in the optically thick limit. The real time half-lives of the exponential fits that result from the 
product of pure ^^Co exponential decay and the f{x) function exponential fits discussed in § 3.2 and put 
on the real time scale with the fiducial core-collapse model are shown in Table 2. All these fits (when 
provided with appropriate constant coefficients) are good fits to the N^^(t) deposition curve in this case. 

Actual large-mass and/or slowly-expanding core-collapse supernovae may never get much beyond the 
early or inflection point 1st order exponential phases of A^]^Jj(t)-like RDE deposition. The reason is that these 
phases could be so long that long-lived radioactive species may become important or the quasi-steady state 
period may end before the phases end. Supernova SN 1987A, for example, which we discuss in § 5 did not get 
past (or much past anyway) the inflection point 1st order phase before both long-lived radioactive species and 
time-dependent effects in the energy processing from decay to UVOIR emission became important. On the 
other hand, small mass and/or quickly expanding core-collapse supernovae (probably most SNe lb and Ic) 
reach lower optical depth sooner and may get well past the first two 1st order exponential phases before 
A''j^i(t)-like RDE deposition ends. In § 6 we consider SN Ic 1998bw as a possible example of this case. 

4. SN la 1992A 

As an example of how far simple RDE deposition calculations can take one in understanding SN la light 
curves we present a comparison in Figures 4a (days — 250) and 4b (days — 1000) between two normalized 
RDE deposition curves calculated using the LS procedure (see § 3.3) for model W7 and the V light curve of 
normal SN la 1992A. The first of the deposition curves (solid line) is a full W7 deposition curve calculated 
with deposition from all the significant radioactive species in model W7 and with X-ray deposition included. 
In Table 3 we summarize the significant radioactive decays for model W7. The other deposition curve 
(dotted line) is just the limited W7 deposition curve presented in § 3.3. This curve is, for our present 
interest, virtually the same as the N^i{t) curve fitted to model W7 (see § 3.3). The full and limited curves 
are in close agreement until very late times (see Fig. 4b). 

At very early times, ''"'Co and ^^Ni contribute significantly to the full curve although this contribution 
is almost invisible on the scale of Figure 4a. (The full curve is normalized by the RDE deposition of the 
limited curve at time zero, and so rises above 1 until about day 1.3.) At time zero, ^^Co and ^^Ni contribute 
about 20 % of the deposition. Their rapid decay rate is the reason for their significant contribution despite 
the fact that they have only trace abundance. But, of course, their rapid decay rate causes their contribution 
to become negligible quickly: down to to ~ 6% by day 3 and ~ 1 % by day 8. After the ^^Co and ^^Ni 
abundances have become exponentially small, ^^Ni and ^®Co, and then '''^Co alone overwhelmingly dominate 
the full curve because ^^Ni is produced in bulk in SNe la unlike the other radioactive species. In model W7 
the initial ^^Ni abundance is 0.58 Mq and this is a typical value calculated for SNe la. The long-lived 
radioactive species, however, become more important as time passes. By day 400 the long-lived radioactive 
species are contributing about 3% and the ^^Co X-rays about 1 % to the full curve RDE deposition. At 
about day 930 after explosion the long-lived radioactive species become dominant: the ^^Co X-rays still 
supplying only about 1 % of the deposition. 

At day 1000 the breakdown of the full curve's deposition is as follows. Photons contribute 29.2 % and PE 
particle kinetic energy 70.8 %. The X-rays contribute 25.9 % and the 7-rays only 3.3 %. The ^®Co, ^''Co, and 
^^Fe contribute 38.9 %, 56.3 %, and 4.6 %, respectively. The ''^Ti and its short-lived daughter ^''Sc contribute 
only 0.16%. Because of ^'^Ti's long half-life, however, "'^Ti/^'^Sc will dominate the RDE deposition by about 
day 5835. 

It should be emphasized that the predictions of radioactive species abundance by model W7 are uncertain 
on the grounds of both nuclear reaction rates and the explosive nucleosynthesis. In the case of ^®Ni, the 
uncertainty is probably small since explosive production of order 0.6 Mq of ^^Ni in normal SNe la is well 
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supported by a variety of empirical evidence. The abundances of trace radioactive species do not have nearly 
such strong support from observations and are much less certain. 

Ideally one would like to compare the RDE deposition curves to SN la UVOIR bolometric light curves 
for the quasi-steady state phase which begins at ~ 60 days after explosion. So far, however, no adequate 
late-time UVOIR bolometric light curves for SNe la have been extracted. The V light curve, however, 
appears to be a reasonable ersatz. It has been suggested that the bolometric correction to the V light curve 
from about 100 days and up to maybe 600 days after explosion is <,0.1 0.2 magnitudes and is relatively 
constant (e.g., Cappellaro et al. 1997). Therefore, as an example, we have fit the V light curve of normal 
SN la 1992A (Cappellaro et al. 1997) to our deposition curves. The V light curve for SN 1992A is one of 
the best observed for SNe la and has the latest observation of a SN la ever: an HST datum from 926 days 
after the B maximum taken by the SINS team (e.g., Kirshner et al. 1993) and reduced by Cappellaro et al. 
(1997). 

To fit the SN 1992A V light curve we have assumed a rise time to maximum light (taken to be 
the B maximum by convention) of 18 days. This rise time is suggested by observational evidence. The 
fairly normal SN la 1990N was discovered at a very early phase 17.5 ± 1 days before maximum light 
(Leibundgut et al. 1991a). And there is empirically estimated rise time to maximum light of 17.6 ±0.5 days 
for SN 1994D (another fairly normal SN la) (Vacca & Leibundgut 1996). Additionally, observations of 
cosmologically remote SNe la (which are usually discovered well before maximum light) suggest a typical 
SN la rise time of ^ 17 18 days (Nugent et al. 1999). Variation from 18 days by a few days are expected 
for some SNe la, but this will not much affect a fit of late-time light curves to theoretical predictions. The 
vertical level of the V light curve was determined by fitting the V light curve by eye to the deposition curves 
in the range from about day 300 to day 400. 

The first thing to notice in Figures 4a and 4b is that the fit from day 51 on is rather good. (The period 
before about day 60 [i.e., before the quasi-steady phase] is that of significant UVOIR diffusion time, and so 
the V light curve in that period is not expected to trace the RDE deposition directly.) The second is, of 
course, the discrepancies in the post-day-51 period. The V data from the days 51 — 109 after the explosion 
are systematically too high by of order 0.1 dex. The V data from days 207 — 451 have a dispersion about 
the RDE deposition curves of order 0.05 dex. The V data point from day 944 (i.e., from 926 days after 
maximum light) is in a phase in which the two RDE deposition curves have diverged. With its estimated 
uncertainty of 0.12 dex, this point is somewhat inconsistent with both of the two deposition curves. 

Reasons for the discrepancies are easily found. First, the V light curve may well not be a sufficiently 
adequate tracer of the UVOIR bolometric light curve. This is probably at least partially the cause of the 
systematic difference between the days 51 — 109 and days 207 — 451 behavior. On the other hand we can 
improve the overall fit by replacing model W7 with fiducial time to = 40.74 days with a model with fiducial 
time to = 47 days: the days 51-109 region is better fit; the days 207-451 region is fit only slightly worse. 
Perhaps, a model with a to value significantly larger than the to value of model W7 is needed. 

The dispersion of the day 207 — 451 data is consistent with the probable uncertainty in these data of 
< 0.1 dex (Cappellaro et al. 1997). 

That the day 944 V datum is close to the RDE deposition curves at all is remarkable. The prediction 
for SNe la is that the infrared catastrophe should set in ^ 400-600 days after explosion (e.g., Axelrod 1980, 
p. 70, 133; Fransson et al. 1996; see also § 2). This means that a V light curve fitted to the RDE deposition 
curve in the early quasi-steady state phase is expected to fall well below the RDE deposition curve later on 
when most of the UVOIR bolometric emission shifts to the infrared. Fransson et al. (1996) suggested that 
perhaps clumping in the ejecta or more optical emission from recombination cascades might prevent the 
infrared catastrophe by keeping strong emission in the optical. In addition to the infrared catastrophe, the 
quasi-steady state should breakdown after about day 600 (e.g., Axelrod 1980, p. 48). A third factor is that 
PE particle escape could reduce the UVOIR bolometric light curve below the complete PE particle trapping 
prediction we have made (see § 2). 

Of course, not too much weight should be placed on a single isolated V datum even if it is accurate. 
However, there is also one very late AB band magnitude (AB — 0.2) for SN la 1972E from day 732 (again 
taking the rise time to maximum as 18 days) that is consistent within uncertainty with quasi-exponential 
decline from the day 98-434 period (Kirshner & Oke 1975). (Note the B light curve is not considered as 
good a tracer of the UVOIR bolometric light curve as the V light curve, but V data beyond day 434 is not 
available for SN 1972E.) Somehow a combination of effects seem to keep the very late-time SN la light curves 
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(out to about day 950) quasi-exponential or nearly so with respect to early late-time behavior. 

We have only considered one SN la here as an example. There certainly are variations between late- 
time V light curves from different SNe la (e.g., Cappellaro et al. 1997). However, there are certain average 
SN la observational results that can be compared to our RDE deposition curves. Leibundgut (1988) and 
Leibundgut ct al. (1991b) presented template light curves for SNe la that extend into in the early quasi- 
steady state phase. The half-life of an exponential fitted to the day 80-120 period of their V template is 29.4 
days. The mean SN la light V light curve of Doggett & Branch (1985) drawn from older data has a half-life 
of 52 days for the period of about 160-320 days after explosion. These half-lives are in fair agreement with 
the half-lives of 29.23 days and 54.95 days for exponential fits to the RDE deposition for the x-rangcs [2, 3] 
(approximately days 80 120) and [4,8] (approximately days 160-320), respectively, for our fiducial SN la 
model (see Table 2). The corresponding half-lives for model W7 N^^{t) curve (which is very close agreement 
with our fiducial model curve) (see § 3.3) are 29.16 days and 54.88 days, respectively. The good agreement 
between the mean observational half-lives and the half-lives of the exponential fits to our fiducial SN la 
model and model W7 N^i{t) RDE deposition curves and the fairly good fit of the model W7 RDE deposition 
curves to the SN 1992A V light curve in Figures 4a and 4b suggest that model W7 and similar models are 
adequate average models for SNe la in respect to RDE deposition. 

Cappellaro et al. (1997) have done a more extensive, detailed analysis of late-time SN la V light curves 
and have included positron transport in their calculations. Probably their most notable conclusion was that 
even normal SNe la could come from a range of masses. The greater masses are needed for the slowest 
declining SNe la: i.e., those whose V light curves approach the late 1st order exponential phase most slowly. 
That their slowest declining SNe la, by one measure at least (viz., the width of nebular phase iron-peak 
element emission lines), seem to have the lowest central concentration of initial ^^Ni makes the need for mass 
variation seem all the stronger. Recall that lower concentration of ^^Ni (implying smaller q and to) aids 
7-ray escape, and so hastens the approach to the late 1st order exponential phase (see § 3.3). 

Cappellaro et al. (1997), however, did not perform a detailed NLTE treatment of the conversion of 
RDE into UVOIR emission. Such treatments have been done (e.g., Axelrod 1980; Ruiz-Lapuente et al. 1995; 
Ruiz-Lapuente 1997; Liu et al. 1997; Liu, Jeffery, & Schultz 1998), but limitations, particularly in atomic 
data input, have limited the conclusions. More detailed, less limited NLTE treatments and better late-time 
observations of SNe la would greatly aid in determining the SN la nature. 

5. SN II 1987A 

Many detailed studies of the late-time emission of SN II 1987A have been done and these have made 
considerable progress in understanding this object (e.g., Woosley 1988; Bouchet et al. 1991; Suntzeff et al. 
1991, 1992; Li, McCray, & Sunyaev 1993; Fransson & Kozma 1993; Nomoto et al. 1994a; Pransson et al. 
1996; Suntzeff 1998, and references therein). Here we simply wish to show that the N^^{t) RDE deposition 
curve for our fiducial core-collapse model (see § 3.4) with fiducial time to = 563.969041 days adequately 
accounts for the early late-time (i.e., the quasi-steady state) UVOIR bolometric light curve of SN 1987A. 

We plot the fiducial curve in Figures 5a and 5b along with the UVOIR bolometric light curve of SN 1987A 
(Bouchet et al. 1991; Suntzeff et al. 1991). We have vertically adjusted the SN 1987A curve to fit the fiducial 
curve in day 130-200 period. The explosion epoch of SN 1987A is, of course, solidly fixed by the SN 1987A 
neutrino burst (Bionta et al. 1987; Hirata et al. 1987). For the day 134-432 period we have only plotted the 
very high accuracy bolometric data derived from spectrophotometry (Bouchet ct al. 1991). 

As can be seen from Figure 5a, SN 1987A very suddenly entered the quasi-steady state phase where 
RDE deposition equals the UVOIR bolometric light curve. The transition of the SN 1987A curve from steep 
decline to slow decline took place over the day 120-128 period (Suntzeff & Bouchet 1991). The early slow 
decline phase of the SN 1987A curve is in fact the manifestation of the early 1st order exponential phase 
of RDE deposition that we discussed in § 3.2. This can be seen by noting that the deposition curve in 
Figure 5a only slowly diverges from the complete trapping curve and only well after the "''^Ni deposition 
phase has ended. Bouchet et al. (1991) give the SN 1987A curve half-life for the day 134-300 period to be 
76.0 ± 0.2 days. For the fiducial model, the day 134-300 period corresponds in reduced time to the x-range 
[0.238, 0.532]. The half-life for an exponential fit to the fiducial ciu'vo for the a;-range [0.23, 0.53] is 75.91 days 
(see Table 2). The observed and fiducial model half-lives are in good agreement. 
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The fiducial curve fits the SN 1987A curve very well from the; beginning of the quasi-steady state phase 
until about day 700. Of course, we chose the parameter values for the fiducial model in order to get this 
good fit (see § 3.4). In reality, however, we have only shown that the N^^it) function with a well chosen 
to value gives a good fit. Thus, the underlying picture of "''^Co RDE deposition in a medium of increasing 
transparency is certainly adequate. However, our parameterized core-collapse model may not itself be very 
adequate to account for SN 1987A. But even if it were exactly right for SN 1987A, we would still not be able 
to determine the right parameter values without more constraints than the SN 1987A curve. Model lOH 
suggests some other constraints, but this is not sufficient to find a unique fit since among other things 
model lOH does not exactly correspond to our parameterized core-collapse supernova model (see § 3.4). 

After ~ 800 days after explosion, the UVOIR bolometric light curve of SN 1987A starts declining 
distinctly less steeply. (This change is not obvious without showing post-day-1000 data.) The change is 
adequately accounted for by ^^Co RDE deposition becoming important (e.g., Suntzeff et al. 1992) and 
the ending of quasi-steady state of the ejecta by the ionization freeze-out (e.g., Fransson & Kozma 1993; 
Fransson et al. 1996). These effects are, of course, not included in the N^i{t) function. 

6. SN Ic 1998bw 

Supernova SN 1998bw is a very interesting object because it is a possible cause of the 7-ray burst 
GRB980425 which occurred on 1998 April 25.91 (Galama et al. 1998; Iwamoto et al. 1998). It was also 
remarkable just as a SN Ic in being unusually bright and having very fast outer ejecta (Iwamoto et al. 1998). 
And the supernova showed net intrinsic polarization (Kay et al. 1988) suggesting large asymmetry (Hoflich, 
Wheeler, & Wang 1998). Intrinsic polarization, however, is not it seems an unusual feature of core-collapse 
supernovae (Wang, Wheeler, & Hoflich 1998). 

Yet another remarkable aspect of SN 1998bw is that its BVI light curves for ~ 60-185 days after 
explosion (almost certainly entirely in the quasi-steady state phase) are very exponential, but do not have 
the ^^Co half-life (McKenzie & Schaefer 1999). (Here we assume that GRB980425 sets the date of explosion. 
From early-time light curve analysis of SN 1998bw, Iwamoto et al. [1998] find that the explosion epoch agrees 
with GRB980425 to within -1-0.7/— 2 days.) As discussed in § 3.2, there is a continuum of quasi-exponential 
phases between the inflection point 1st order exponential phase and the late 1st order exponential phase as 
t ^ 00. The half- lives determined for exponential fits to the riuasi-cxponential phases will not be the ^^Co 
half-life. At first one might suppose that the SN 1998bw light curves for the day 60-185 period are from 
a quasi-exponential phase. The light curves, however, appear truly exponential, not just quasi-exponential. 
Over the day 60 185 period, the data with uncertainty of <,0.1 magnitudes (^0.04 dex) shows perhaps 
only a slight trace of curvature on a semi-logarithmic plot and that trace has been deemed insignificant 
(McKenzie & Schaefer 1999, Fig. 1). 

It may be in the SN 1998bw case that the light curves are showing the signature of the inflection point 
1st order exponential phase of RDE deposition that occurs near t = to (see § 3.2). Militating against this 
idea is the fact that the light curve half-lives are not all the same. The B, V, and I half-lives are 53.4 ± 0.8, 
40.9 ± 0.7, and 41.6 ± 0.7 days, respectively. If V and / dominated the emission, then one could tentatively 
conclude that they trace a UVOIR bolometric light curve with a half-life of about 41 days. But spectra 
through day 136 show that the flux in the B band is comparable to that in the VI bands (Patat et al. 
1999). The sum of different exponentials is not exactly exponential. Nevertheless it is difficult to believe 
that one can get close to exponentials in three bands, two of them with almost the same half-life, without 
the driving RDE deposition being nearly an exponential. We wiU assume that the UVOIR bolometric light 
curve for the day 60-185 period is nearly an exponential with the UVOIR bolometric light curve half-life 
44 days estimated by McKenzie & Schaefer (1999). 

Figure 6a shows a schematic SN 1998bw UVOIR bolometric light curve for the day 60 185 period 
with half-life 44 days, a fitted N^^it) RDE deposition curve, and the complete trapping curve. The N^i{t) 
curve has been fitted to the schematic curve by the following procedure. The vertical level of the schematic 
SN 1998bw curve was determined by a least-squares fit to the N^[{t) curve. The N^[{t) curve was then 
varied by varying to (which is N^^^tys only free parameter when rch,o is set to 1 as it always is) in order 
to minimize the maximum deviation between the N^^{t) curve and the schematic SN 1998bw curve. The 
final fitted N^^it) curve has a fiducial time to = 134.42 days (for Tch,o = !)• The maximum deviation of 
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fit is 0.024 dcx (0.06 magnitudes). McKcnzic & Schaefer (1999) find no evidence for systematic deviations 
from exponential behavior in tfie broad band light curves over the day 60-185 period greater than 0.02 dex 
(0.05 magnitudes). If the actual UVOIR bolometric light curve of SN 1998bw has the same closeness to exact 
exponential behavior as the observed broad band light curves, then our fitted N^-^{t) curve has about as much 
deviation from exponential behavior has can be tolerated. The maximum deviations of iVjJj(t) curves with 
fiducial times 90, 120, 150, and 180 days are 0.12 dex, 0.048 dex, 0.061 dex, and 0.11 dex, respectively. None 
of these curves could be tolerated given the supposed exponential behavior of the UVOIR bolometric light 
curve. 

Our fitted to = 134.42 days puts the inflection point of the RDE deposition at day 139.90 (see § 3.2). 
Thus, as we could have anticipated, a best fit to a highly exponential region of UVOIR bolometric light 
curve that does not have the ^^Co half-life can force the fitted region to be in the inflection point 1st order 
exponential region of the N^i{t) function. That the N^i{t) approximation tends to be weakest when time 
is of order to (see § 3.2) is a weakness of our analysis. But a single inflection point 1st order exponential 
region is expected even in a more sophisticated treatment of RDE deposition provided the supernova ejecta 
is fairly smoothly varying and the ^^Ni contribution becomes exponentially small soon enough (see § 3.2). 

In Figure 6b we predict RDE deposition curve and quasi-steady state UVOIR bolometric light curve out 
to day 1000 using N^^it) with to = 134.42 days. The schematic SN 1998bw curve is again displayed. The 
actual RDE deposition curve will probably differ from the prediction due to long-lived radioactive species; 
the UVOIR bolometric light curve due to those species and time-dependent effects (see § 2). The infrared 
catastrophe (see § 2) does not in itself cause deviations from A7'];^;(t)-like behavior in the UVOIR bolometric 
light curve, but it does make it harder to measure the UVOIR bolometric light curve. The predicted RDE 
deposition curve's approach to the complete escape curve is sufficiently slow that non-exponential behavior 
in the UVOIR bolometric light curve would be difficult to detect without high quality data out to several 
hundreds of days after explosion. 

Taking our fitted N^^{t) curve at face value, we can make an estimate of the ejecta mass by inverting 
an expression for tg for a given parameterized model. We will assume the parameterized core-collapse model 
of § 3.4 and invert equation (33) for mass. Recall that this parameterized model assumes that the bulk of 
the ejecta has density going as the inverse square of velocity. Some support for this density profile for outer 
ejecta is provided by Branch (1999) who found that it worked well in a parameterized LTE analysis of the 
photospheric epoch spectra of SN 1998bw. 

We now need to estimate the model parameters. From Iwamoto et al. (1998) it is known that 
SN 1998bw ejecta extend out to velocities > 28,000kms~^. The analysis of Branch (1999) suggests ejecta 
velocities >, 60, 000 kms~^ and that the inverse-square density profile extends even that far. We will be 
more conservative and assume that the vast bulk of the ejecta is confined to an inner high-mass core as 
in model lOH upon which we based our parameterized core-collapse model. The estimated half-width of 
one emission line of the [O I] AA6300,6364 emission above the estimated pseudo-continuum of the day 136 
nebular spectrum of SN 1998bw (Patat et al. 1999) corresponds to a velocity of order 20000 km s~^. We 
assume 20000 km s~^ is the outer velocity Vh of the high-mass core. Given that the analysis of Branch (1999) 
and that the pseudo-continuum is ill-defined, 20000 km s~^ could be a factor of 2 or 3 too small. For the 
inner velocity of the liigli-mass core, we have no strong evidence. But if there was a large hollow at the 
center of the ejecta the unblended emission lines would tend to be fiat-topped (e.g., Mihalas 1978, p. 477; 
Jeffery & Branch 1990, p. 190). The [O I] AA6300,6364 emission in the day 136 spectrum has a very sharp 
peak (which is probably mostly due to the strongest line), and we estimate that a low-mass, low-density 
center effectively setting Wa >, 2000 kms~-^ is unlikely. We will just assume = 700 km s^^ as we did for 
our fiducial core-collapse model. For the concentration factor q we choose 1 implying that the "''^'Ni was 
concentrated in the innermost regions of the high-mass core. Note that Iwamoto et al. (1998) concluded 
that large-scale mixing of the ^^Ni was needed to fit the fast rise of the estimated pre-maximimi SN 1998bw 
UVOIR bolometric light curve. The real q could be significantly less than 1. For k we choose 0.028 cm^g^^. 
This an appropriate value for all-metal ejecta in the transition phase between optically thick and thin limits 
as determined from Table 1 of Jeffery (1998a) or Table 2 of Jeffery (1998b). (SNe Ic are believed to have no 
significant hydrogen or helium.) In extracting our to value from the fit to the schematic SN 1998bw curve 
we have already set Tch,o = 1 as we always do. 
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The expression for mass for our parameterized core-collapse model is 



M =15 Mq X 



tn 



563.969041 days 

1 - {Va/Vb) 



0.035 cm^g- 



1 - {v./vb) (1/g) 



( ^2 ) 

V700 knis-iy 



5000 kms- 



(34) 



With the chosen parameter values, we obtain a mass of 4.26 Mq. Even given that our parameterized core- 
collapse model is appropriate for SN 1998bw, this value could be factor of few too small: could be 2 or 3 
times what wc assiimc and q coiild significantly less than 1. The optimum value could also be different 
from what we have assumed by a factor of 2 or more either up or down. But since our parameterized 
core-collapse model may not be appropriate, our mass estimate may be only order of magnitude in any case. 

Other mass determinations or lower limits on mass for SN 1998bw have been given. For example, 
Iwamoto et al. (1998) find that a 13.8 carbon-oxygen star explosion model can reproduce fairly well 
estimated early-time SN 1998bw UVOIR bolometric light curve. And Branch (1999) from his photospheric 
epoch parameterized LTE analysis estimated that SN 1998bw had ^ 6 Mq above 7000 km s"-'^. There is still, 
however, considerable debate about the SN 1998bw mass and structure. Given an adequate parameterized 
structural model for the supernova, the to parameter obtained by fitting the UVOIR bolometric light curve 
from the quasi-steady state phase would help to constrain the supernova's mass and/or other parameters. 
Since SN 1998bw may be highly asymmetric, the determination of an adequate parameterized structural 
model may be diflacult. 



7. CONCLUSIONS 

We have given a presentation of a simple, approximate, analytic treatment of RDE deposition 
in supernovae from the decay chain ^^Ni — >^^Co — >^^Fe. The treatment provides a straightforward 
understanding of the exponential/quasi-exponential behavior of the UVOIR bolometric luminosity and a 
partial understanding of the exponential/quasi-exponential behavior of the broad band light curves. The 
treatment reduces to using the normalized N^i{t) deposition function (see § 3.2) as an analysis tool. (The 
absolute deposition is determined by specifying the initial ^^Ni mass or fitting absolute UVOIR supernova 
luminosity.) The time evolution of N-^^(t) is determined by three time scales: the half-lives of ^^Ni and 
^^Co, and a fiducial time parameter to that governs the 7-ray optical depth behavior of a supernova. The 
to parameter can be extracted from a structural supernova model, and we have shown examples of how 
this is done in §§ 3.3 and 3.4. It can also be obtained by fitting the RDE deposition curve from a more 
detailed treatment of deposition as we did in § 3.3 or by fitting to an observed UVOIR bolometric light curve 
from the quasi-steady state phase of a supernova as we did in §§ 5 and 6. A Iq parameter extracted from 
observations can provide a constraint on the important physical parameters of a supernova. Effective use of 
this constraint, however, requires having an adequate parameterized structural supernova model. 

The N-^^{t) function is used to analyze the preliminary UVOIR bolometric light curve of SN 1998bw 
(the possible cause of GRB980425) (§ 6). The SN 1998bw fiducial time to is found to be 134.42 days and a 
prediction is made for the evolution of the SN 1998bw RDE deposition curve and quasi-steady state UVOIR 
bolometric light curve out to day 1000 after the explosion. A crude estimate (perhaps a factor of a few too 
small) of the SN 1998bw mass obtained from a parameterized core-collapse model and to — 134.42 days is 
4.26 Mq. As further examples of the simple analytic treatment, the RDE deposition and luminosity evolution 
of SN la 1992A and SN II 1987A have also been examined (see §§ 4 and 5). 

The simple analytic treatment of RDE deposition has actually existed at least since Colgate et al. 
(1980a, b), but has not hitherto been given a detailed or general presentation as far as we know. The main 
value of this paper is the explicit, detailed, general presentation of this analytic treatment. 
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APPENDIX A 

THE EXPONENTIAL DENSITY PROFILE MODEL 
FOR HOMOLOGOUSLY EXPANDING SUPERNOVAE 

As discussed in § 3.3, hydrodynamic explosion models for SNe la often have density profiles in the 
homologous expansion epoch that are close to exponentials (i.e., inverse exponentials) as functions of radial 
velocity. In this appendix wc present some useful analytic results for exponential density profile models for 
homologously expanding supernovae. 

First, we note the following general integral solution that is useful in developing the analytic results: 



Uz) = / dz' {z'T exp i-z') = exp {-z) ^ ^z'^ , (Al) 

"^^ fe=0 



where n > is an integer. 

Next recall that the radius r of a mass element in homologous expansion after the initial radii have 
become insignificant is given by 

r = vt, (A2) 

where v is the mass element's radial velocity and t is the time since explosion. Recall also that the element's 
density at any velocity declines as t~^. We will use radial velocity as comoving radial coordinate and define 
a dimensionless radial coordinate z by 

z = v/ve , (A3) 

where Ve is the e-folding velocity of an exponential model's density profile. The expression for the density 
profile can then be written 



p{v,t) = pce,o j exp{-z) , (A4) 

where Pcc.o is the central density at fiducial time to- 

Using the equations (A1)-(A4), the expression for mass exterior to radius z for an exponential model is 

M{z) = 4wt^J dvv^p{v,t) = Mexp{-z) ^1 + z+^z'^^ , (A5) 

where M = Snpcefi {veto)^ is total mass. Using the same equations, the expression for kinetic energy exterior 
to radius z is 



E{z) = Ant^ / dv —p{v, t) 



6Mt;2exp(-^)( 1 + ^+1^2 ^1^3 ^1^4) ^ (^6) 

where E = 487rpce,o (we^o)'' = &Mv1 is total kinetic energy. 

It is often useful to have expressions for central density pco.o, total kinetic energy E, and the ratio of 
total kinetic energy to total mass E/M in terms of fiducial parameter values. Since model W7 (Nomoto, 
Thielemann, & Yokoi 1984; Thielemann, Nomoto, & Yokoi 1986) is a widely-used standard SN la model 
whose density closely approximates an exponential, we will use it as a basis for the fiducial values. Model W7 
is a Chandrasekhar mass model, and so has total mass — 1.38 Mq. The equivalent-exponential model Vg 
for model W7 in round numbers is 2700 km s^^. (Wc explain equivalent-exponential model below.) For a 
fiducial time to, we choose 1 day for convenience in simple calculations. (Note this fiducial time is not for the 
same purpose as the fiducial time we use in the main text for optical depth evolution.) Using M = 1.38 M©, 
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2700 km s ^, and to = I day as fiducial values, the desired expressions are 

M 



Pce,0 



0.860327 X 10-**gcm-^ x 



M 



, 1.38 Mq 
/2700kins-i\^ /Iday^ ^ 



to 

£; = 6M.- 1.20064foex{^)(^^^5^) , (A8) 



(A7) 



and 

E 
M 



6..^ = 4.37400 X 10- crgsg-x(^^^^)^ 

0.870032foeM-x(^^^^)\ (A9) 



where a foe (for ten to the fifty-one ergs) is a standard supernova energy unit of 10^^ ergs. (Recall from the 
main text that for numerical consistency wc treat fiducial values and the solar mass unit Mq = 1.9891 x 10^^ g 
[Lide & Frederikse 1994, p. 14-2] as nearly exact numbers.) 

Since SN la hydrodynamic explosion models often have density profiles that are close to exponential, 
there is some interest in specifying exponential models that in some ways arc equivalent to those explosion 
models: viz., equivalent-exponential models. Two parameters are needed to specify an equivalent-exponential 
model. The explosion model's central density and an e-folding velocity drawn from a fit to the explosion 
model's density profile may not be good choices for these parameters: the central density may not be 
representative of the overall explosion model and good criteria for obtaining the fitted e-folding velocity 
need to be specified. We will instead choose the explosion model's total mass M and total kinetic energy E 
as the parameters. Together these parameters should yield an equivalent-exponential model that globally is 
much like the original explosion model: the closer the explosion model is to an exponential model, the better 
the likeness, of course. The expression for the equivalent-exponential model e-folding velocity follows from 
equation (A8) or equation (9): 



= Jl^ = 2700kms- X , ( ^ , ) (1:^) . (AlO) 
^ V6M Y V 1.20064 foey V M J ' 

Given Vg, the central density of the equivalent-exponential model is obtainable from equation (A7). 

We have calculated the equivalent-exponential model Ve for model W7 to be 2724 km s^^. Since the 
calculation depends a bit on how one integrates the kinetic energy over the finite number of zones that make 
up model W7, there is no real reason for using exactly 2724 km s~^ as a fiducial value. Thus, above we 
chose the round number 2700 km s~^. 

It is useful to specify a few other exponential model results. The central atom density riat is given by 

Pcc,0 

na.t = - 



^amuMat 

0.925179 X 10" cm"^ 



(r^) ('J^!^) (l^) . (All) 



/iat/ \1.38MqJ V Ve J \ to 

where mamu = 1.6605402(10) x lO^^'^g (Lide & Frederikse 1994, p. 1-1: uncertainty in the last digits is 
given in the brackets) is the atomic mass unit (amu) and is the mean atomic mass, and where the second 
expression is in terms of fiducial values. The mean atomic mass is defined by 
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where the sum is over all elements ?', Xi is the mass fraction of clement i, and Ai is element i's atomic mass. 
Since the center of model W7 is dominated by stable ^^Fe (even at time zero) and this may be typical of 
SNe la, we chose the fiducial value of /Xat to be 56, the whole number atomic mass of ^^Fe. Note that even 
if ^^Fe (or ''^Ni or ''^Co) do not dominate the center, iron peak elements (with A values fairly close to 56) 
almost certainly do. 

The central free electron density is given by 

Pce.O 



Ue = 



0.925179 X IQi'^cm-^ 

56^/ M \ / 2700km s-i\Vl day V 



^ilj yi.SSMQj V Ve J \ to J ' ^^^^^ 

where /!,, is the mean atomic mass per free electron and where the second expression is again in terms of 
fiducial values. The mean atomic mass per free electron is defined by 

i 

where the sum is over all ions i, is the mass fraction of ion i, Zi is the charge on ion i, and Ai is ion I's 
atomic mass. Note that Jle is not the same as fie, the mean atomic mass per electron defined by equation (10) 
in § 3.1: jjLg accounts for all electrons, not just free electrons. The ionization stage of the center of ejecta 
probably varies strongly as a function of time. But since the central iron-peak elements are likely to be at 
least singly ionized until very late times, it seems most convenient just to choose the singly ionized state of 
an ^ = 56 species as the fiducial ionization. Calculations suggest that the central iron will be mostly singly 
ionized at about day 300 after explosion (e.g., Ruiz-Lapuente et al. 1995; Liu et al. 1998). 
For a constant opacity k, the radial optical depth from radius z to infinity is 

t[z) = KPcefiVet {^J J eJCp{-z) = ^^^2f2 exp(-2) 

= Tcefl(^-j-] exp(-z) , (A15) 



where rce,oi defined by 



kM 

is the radial optical depth to the center at the fiducial time to . The constant opacity 7-ray and free electron 
(i.e., Thomson) radial optical depths from the center to infinity in terms of the fiducial values are 

5017.42 



'ce.O 



M 



0.025 cm2g-i J V1.38Mq 
/'2700kms-iy /I day 



(A17) 



and 



1435.77x('^^^ 



fie J V 1-38 Me 

respectively. For the fiducial value of k for 7-rays we chose 0.025 cm^ g"^ which is a good general value for the 
effective absorption opacity for ^^Co 7-rays for an all-metal medium in the optically thin limit (Swartz et al. 
1995; Jeffery 1998a, b; see also § 3.1). The electron opacity is given by 

(Tp 0.40062033 
Ke = — = z , (A19) 
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where (Te = 0.66524616(18) x 10^24 ^,^^2 ^Lidc & Frcderikse 1994, p. 1-2 : uncertainty in the last digits is given 
in the brackets) is the Thomson cross section. For the fiducial value of /ig we again chose 56 for niceness even 
though the ionization state and composition of ejecta vary widely with velocity location and the ionization 
state with time also. 
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FIGURE LEGENDS 



FIG. la. The N^^{t) deposition curves for a range of to values. The complete trapping curve is a N-^^{t) 
curve with to = co. The complete escape curve is a N-^^{t) curve with to = 0. The complete trapping and 
escape curves are shown on all subsequent deposition curve figures for convenient reference. 
FIG. lb. — The absorption function f{x) (for ^^Co) and exponential fits to its inflection point, the point at 
infinity, and the x-range [2,3]. The fit to the point at infinity is just the constant asymptote which f[x) 
approaches as a; — > oo. 

FIG. 2a. -The normalized deposition curve for model W7 calculated with only ^^Ni and '"'^Co and no X-rays 

(i.e., the limited W7 deposition curve) and the fitted analytic iV^i(^) W7 deposition curve. 

FIG. 2b. — The same as Fig. 2a, except extending to day 500. 

FIG. 3. — The N^^{t) fiducial core-collapse supernova (SN) deposition curve. 

FIG. 4a. — The normalized deposition curve for model W7 calculated with all important radioactive species 
and X-rays included (i.e., the full W7 deposition curve), the limited W7 deposition curve, and the V light 
curve of SN la 1992A. The V light curve has been vertically shifted to fit the deposition curves in roughly 
the day 300-400 period (which is shown in Fig. 4b). The fit is simply determined by eye. 
FIG. 4b. — The same as Fig. 4a, but extending to day 1000. 

FIG. 5a. — The UVOIR bolometric Hght curve of SN II 1987A compared to the N^i{t) fiducial core-collapse 
supernova (SN) deposition curve. The SN 1987A curve has been vertically shifted to fit the deposition curve 
in the day 130-200 period. The fit is simply determined by eye. 
FIG. 5b. — The same as Fig. 5a, but extending to day 1000. 

FIG. 6a.— The schematic SN Ic 1998bw UVOIR bolometric Ught curve (half-life 44 days) and a fitted A^Ni(*) 

deposition curve with fiducial time to = 134.42 days. 

FIG. 6b. — The same as Fig. 6a, except extending to day 1000. 
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TABLES 



TABLE 1 

PARAMETERS FOR THE RADIOACTIVE DECAYS 
OF 56Co AND 56Ni 



Par;uuot(^r 


'■"Co 


'«Xi 


ti/2 (days) 


77.27(3) 


6.077(12) 


te = ii/2/ln(2) (days) 


111.48(4) 


8.767(17) 


Ototai (MeV) 


4.566(2) 


2.135(11) 


Oph+PE (MeV) 


3.74(4) 


1.729(17) 


Qph (MeV) 


3.62(4) 


1.723(17) 


/ph 


0.968(14) 


0.996(14) 


/PE 


3.20(16)-2 


3.99(18)-3 


C (ergss~^ g~^) 


6.70(7)+9 


3.94(4)+10 


B (ergss~^g~^) 


7.27(7)+9 





NOTE. — The values have been taken or derived from LBL, Huo 1992, and Browne & Firestone 1986. 
We have put the uncertainties in the last digits of the parameters in brackets and have written x 10^*^ as 
±k. 

The parameters are defined as follows: ti/2 is half-life, te is e-folding time, Qtotai is the total energy 
(including neutrino energy) per decay, Qph+PE is the mean photon phis PE (positron and dccay-cjcctcd 
atomic electron) kinetic energy per decay, Qph is the mean photon energy per decay, /ph is the fraction of 
<3ph+PE in photon energy, /pe is the fraction of Qph+PE in PE kinetic energy, and C and B are energy 
generation coefficients specified in the text (see § 3.1). The 7-ray energy from positron annihilation is 
included in Qph+PE and Qph. We assume that the neutrinos simply escape the supernova ejecta and make 
no contribution to the RDE deposition. 
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TABLE 2 
EXPONENTIAL FITS 
TO THE ABSORPTION FUNCTION f{x) 



Fit at x/over [xajXb] 

Max. Error (dex) in [a;a,a;b] 


-^coef 


Xl/2 


.fid, la 
''1/2 

(days) 


.fid,CC 
''1/2 

(days) 


x = 

0.001 in [0, 0.40] 
0.04 in [0, 0.64] 


1 


oo 


77.27 


77.27 


[0.23,0.53] 
0.0033 in [0.23,0.53] 


1.027959 


7.659517 


61.98 


75.91 


X = a;infi ~ 1.040765 
0.004 in [0.8,1.35] 
0.04 in [0.56,1.82] 


1.950777 


0.625347 


19.21 


63.38 


[0.54, 1.98] about a;infl (n) 
0.014 in [0.54,1.98] 


1.694608 


0.705462 


21.01 


64.70 


[0.42, 2.46] about XinS (m) 
0.028 in [0.42,2.46] 


1.540518 


0.770151 


22.38 


65.60 


[0.33, 2.82] about Xinfi (b) 
0.04 in [0.33,2.82] 


1.441097 


0.822176 


23.43 


66.23 


[2,3] 
0.01 in [2, 3] 


0.804766 


1.149675 


29.23 


69.04 


[4,8] 
0.037 in [4, 8] 


0.151079 


4.653005 


54.95 


75.06 


a; = 00 

0.098007 in [2a;ti.,cxD] 


0.0320 


CXD 


77.27 


77.27 



NOTE. — The exponential function Xcoof exp {—x/xe) (where Xe is the e-folding reduced time) has been 
fitted to the absorption function f{x). The fits are either to /(x)'s value and its logarithmic slope at a given 
point X or are fits to f(x) over a given x-range [xa,a;b]. The point fits are to the asymptotic behavior of 
f{x) at zero and infinity, and to f{x) at its inflection point Xinfl. A fit to an x-range was chosen so as to 
reproduce /(x) close to optimally over that range. 

The first column gives the point or x-range of the fit in the first line and in subsequent lines the 
maximum error of the fit (or the maximum deviation of /(x) from the fitted exponential) in specified ranges. 
The second column gives the coefficient Kcoef of the fitted exponential. The third column gives the reduced 
time half-life of the exponential: Xi/2 = Xeln(2). The fourth and fifth columns give the half-lives (in real 
time) of the product of the ^®Co decay exponential and the fitted exponential for the fiducial SN la model 
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and fiducial core-collapse (CC) supernova model posited in §§ 3.3 and 3.4, respectively. 
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TABLE 3 

SUMMARY OF RADIOACTIVE DECAYS IMPORTANT IN MODEL W7 



Radioactive species and initial mass 



Decay 



t 



1/2 



Qtotal 

(MeV) 



Qph+PE 

(MeV) 



Qph 

(MeV) 



0) 



44rpj _^44sj. 
44Sc ^44ca 



"^"Sc 1.8 X 10-9 M« 



63(3) years 
0.1636(3) days 



0.2675(19) 0.1493(15) 0.1384(14) 

3.6533(19) 2.73(23) 2.13(23) 



s^Co 4.5 X 10-3 Mq; 



55Co ^55Fg 
55pg ^55]y[j^ 



55 



Fe 1.4 X 10-3 Mq 
0.7304(13) days 
2.73(3) years 



3.4513(4) 
0.23138(10) 



2.43(4) 
0.0054(3) 



2.00(4) 
0.00163(5) 



56Ni 0.58 M©; ^'^Co 6.1xlO-'^Afo 
56Ni ^56cq 6.077(12) days 

56Co ^56pg 77.27(3) days 



2.135(11) 
4.566(2) 



1.729(17) 
3.74(4) 



1.723(17) 
3.62(4) 



57 



Ni 2.15 X 10-2 Mq; 

57Co ^57Fe 



5^Co 8.1 X 10-^ Mq 
1.4833(25) days 3.264(3) 
271.79(9) days 0.8360(4) 



2.07(3) 
0.1429(8) 



1.92(3) 
0.1253(6) 



NOTE. — The masses of the radioactive species are from epoch just after the model W7 explosion (i.e., 
effectively time zero) (Thielemann et al. 1986; Nomoto et al. 1994b). The nuclear data have been taken or 
derived from LBL, Huo 1992, Bhat 1992, and Browne & Firestone 1986. We have put the uncertainties in 

the last digits of the parameters in brackets. 

See the note to Table 1 for the definitions of the parameters. 
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